1 



Chapter 1 

Black Hole Perturbation 

Yasushi Mino 1,2 , Misao Sasaki 1 , Masaru Shibata 1 , 
Hideyuki Tagoshi 3 and Takahiro Tanaka 1 

1 Department of Earth and Space Science, Osaka University, Toyonaka 560, Japan 
' 2 Department of Physics, Kyoto University, Kyoto 606, Japan 

Q\ • ^National Astronomical Observatory, Mitaka, Tokyo 181, Japan 

0\ . 

Abstract 

' In this chapter, we present analytic calculations of gravitational waves from a 

Q ■ particle orbiting a black hole. We first review the Teukolsky formalism for dealing 
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with the gravitational perturbation of a black hole. Then we develop a systematic 
method to calculate higher order post-Newtonian corrections to the gravitational 
waves emitted by an orbiting particle. As applications of this method, we consider 
orbits that are nearly circular, including exactly circular ones, slightly eccentric ones 
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f — ' and slightly inclined orbits off the equatorial plane of a Kerr black hole and give 

the energy flux and angular momentum flux formulas at infinity with higher order 
post-Newtonian corrections. Using a different method that makes use of an analytic 
series representation of the solution of the Teukolsky equation, we also give a post- 
Newtonian expanded formula for the energy flux absorbed by a Kerr black hole for 
a circular orbit. 
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§1. Introduction 



In this chapter, we review recent progress in the analytic calculations of gravita- 
tional waves from a particle orbiting a black hole using a systematic post-Newtonian 
expansion method. There has been substantial activities in this field recently and 
there is a diversity of literature. Here we are mostly concerned with the actual cal- 
culations of the gravitational waves from an orbiting particle and we intend to make 
this chapter as self-contained as possible. We do not, however, discuss much about 
implications of the results to actual astrophysical situations. 

In the black hole perturbation approach, one considers gravitational waves from 
a particle of mass fi orbiting a black hole of mass M assuming [i <C M. Although 
this method is restricted to the case when \i <C M, one can calculate very high 
order post-Newtonian corrections to gravitational waves using a relatively simple 
algorithm in contrast with the standard post-Newtonian analysis. This is because the 
fully relativistic effect of the spacetime curvature is naturally taken into account in 
the basic perturbation equation. We can also calculate numerically the gravitational 
waves without assuming the slow motion of its source. Then, we can easily investigate 
the convergence of the post-Newtonian expansion by comparing the result of the post- 
Newtonian approximation with the fully relativistic one. In this sense, the black hole 
perturbation method gives a very important test of the post-Newtonian expansion. 
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Further, since the effect of the spacetime curvature is naturally taken into account, 
we can easily investigate interesting relativistic effects such as tails of gravitational 
waves. 

We consider the post-Newtonian wave forms and luminosity which are expanded 
by v/c, where v is the order of the orbital velocity. The lowest order of the gravi- 
tational waves are given by the Newtonian quadrupole formula. We call the post- 
Newtonian formulas for the wave forms and luminosity which contain terms up to 
0((v/c) n ) beyond the Newtonian quadrapole formula as the n/2PN formulas. 

Let us first briefly give a historical review. The gravitational perturbation equa- 
tion of a black hole using the Newman- Penrose formalism© was derived by Bardeen 
and Presstf for the Schwarzschild black hole, and by TeukolskycP for the Kerr black 
hole. By using these equations, many numerical calculations of gravitational waves 
induced by the presence of a test particle have been done. We do not list up all 
such works. Here, we only referto three articles; BreuercP, ChandrasekharEP, and 
Nakamura, Oohara, and KojimacP. 

On the other hand, analytic calculations of gravitational waves produced by the 
motion of a test particle have not been developed very much until recently. This 
direction of research was first done by Gal'tsov, Matiukhin and PetukhovlZP in which 
they considered a case when a particle moves a slightly eccentric orbit around a 
Schwarzschild black hole, and calculated the gravitational waves up to 1PN order. 
Then, Poissonlt* considered a case of circular orbit around a Schwarzschild black hole 
and calculated the wave forms and luminosity to 1.5PN order at which the tail effect 
appears. Cutler, Finn, Poisson and SussmancP also worked on the same problem 
numerically by using the least square fitting method, and obtained a formula for the 
luminosity to 2.5PN order. Subsequently, a highly accurate numerical calculation 
was carried out by Tagoshi and Nakamura They obtained the formulas for the 
luminosity to 4PN order numerically by using the least square fitting method. They 
found the logv terms in the luminosity formula at 3PN and 4PN orders. They 
showed that, although the convergence of the post-Newtonian expansion is slow, the 
luminosity formula which is accurate to 3.5PN order will be good enough to represent 
the orbital phase evolution of coalescing compact binaries accurately. After that, 
Sasaki 0^ found an analytic method and obtained the formulas which are needed 
to calculate the gravitational waves to 4PN order. Then, Tagoshi and Sasakill3 
obtained the gravitational wave forms and luminosity to 4PN order analytically, and 
confirmed the results of Tagoshi and Nakamura^ These calculations were extended 
to 5.5PN order by Tanaka, Tagoshi, and SasakiEJ). 

In the case of orbits around a Kerr black hole, Poisson calculated the 1.5PN 
order corrections to the wave forms and luminosity due to the rotation of the black 
hole and showed that the result agrees with the standard post-Newtonian effect due 
to spin-orbit couplingta). Then, Shibata, Sasaki, Tagoshi and Tanakali^ calculated 
the luminosity to 2.5PN order. They calculated the luminosity from a particle in 
circular orbit with small inclination from the equatorial plane. They used the Sasaki- 
Nakamura equation as well as the Teukolsky equation., This analysis was extended to 
4PN order by Tagoshi, Shibata, Tanaka and SasakiliU) in which the orbit of the test 
particle was restricted to circular ones on the equatorial plane. The analysis in the 
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case of slightly eccentric orbit on the equatorial plane was also done by Tagoshi IIZP 
to 2.5PN order. 

Tanaka, Mino, Sasaki, and ShibataEJ) considered the case when a spinning par- 
ticle moves a circular orbits near the equatorial plane around a Kerr black hole, and 
derived the luminosity formula to 2.5PN order including the linear order effect of 
the particle's spin. They used the equations of motion of Papapelrou'sliS) and the 
energy momentum tensor of the spinning particle given by DixonE^P. 

The absorption of gravitational waves into the black hole horizon, appearing 
at 4PN order in the Schwarzschild case, was calculated by Poisson and Sasaki in 
the case when a test particle is in a circular orbit §>. The black hole absorption in 
the case of rotating black hole appears at 2.5PN order c$. Recently a new analytic 
method to solve the homogeneous Teukolsky equation was found by Mano, Suzuki, 
and TakasugiEi 1 . Using this method, the black hole absorption in the case of rotating 
black hole was calculated by Tagoshi, Mano, and TakasugiEJ) to 6.5PN order beyond 
the quadrupole formula. 

If gravity is not described by the Einstein theory but by the Brans-Dicke the- 
ory, there will appear scalar type gravitational waves as well as transverse-traceless 
gravitational waves. .-Such scalar type gravitational waves were calculated by Ohashi, 
Tagoshi and Sasaki in a case when a compact star is in a circular orbit on the 
equatorial plane around a Kerr black hole. 

The organization of this chapter is as follows. We review the Teukolsky for- 
malism for the black hole perturbation in section ^ and formulate a post-Newtonian 
expansion method of the Teukolsky equation in section ^. Then we turn to the 
evaluation of gravitational waves by an orbiting particle in the rest of sections. 

First we consider circular orbits. In section we calculate the gravitational wave 
luminosity from a test particle in circular orbit arouncLa Schwarzschild black hole 
to 5.5PN order, based on Tanaka, Tagoshi, and Sasaki tS'. This is the highest post- 
Newtonian order achieved so far. Based on this result, we investigate the convergence 
property of the post-Newtonian expansion in section ||. In section |6|, we consider 
circular orbits on the equatorial plane around a Kerr black hole and calculate the 
luminosity to 4PN order, based on Tagoshi, Shibata, Tanaka, and SasakiE^. We 
find the luminosity contains the terms which describe the effect of not only spin- 
orbit coupling but also the effect of higher multipole moments of the Kerr black 
hole. 

Next we consider slightly noncircular orbits. In section [?], we calculate the 
0(e 2 ) corrections to the 4PN energy and angular momentum flux formulas in the 
case of a slightly eccentric orbit around a Schwarzschild black hole, where e is the 
eccentricity. In section ||, we consider a slightly eccentric orbit on the equatorial 
plane around a Kerr black hole and evaluate the 0(e 2 ) corrections to 2.5PN order, 
based on TagoshiliZP . Then in section ||, we calculate the gravitational waves induced 
by a test particle in circular orbit with small inclination from the equatorial plane 
around a Kerr black hole and evaluate the 2.5PN_energy and angular momentum 
fluxes, based on Shibata, Sasaki, Tagoshi, Tanakali^P. In section 1C, we discuss the 
adiabatic orbital evolution around a Kerr black hole under radiation reaction and 
show that circular orbits will remain circular under adiabatic radiation reaction but 



4 



the stability of circular orbits can only be examined by an explicit evaluation of the 
backreaction force. 



In section |ll], we consider the effect of the spin of a particle. We first give 
a general formalism to treat the gravitational radiation from a spinning particle 
orbiting a Kerr black hole. Then we calculate the 2.5PN luminosity formula with 
the first order corrections of the spin for circular orbits which are slightly inclined 
due to the spin of the particle. 



Finally, in section 12, we review a calculation of the black hole absorption based 
on Tagoshi, Mano, and TakasugiBS*. The black hole absorption effect appears at 
0(v 5 ) relative to the Newtonian quadrupole luminosity for a Kerr black hole, while 
at 0(i> 8 ) for a Schwarzschild black hole. We show the energy absorption rate to 
0(v 8 ) beyond the lowest order for the Kerr case, i.e., 0(v 13 ) or 6.5PN order beyond 
the Newtonian quadrupole luminosity. 

Since many of the calculations encountered in this black hole perturbation ap- 
proach are lengthy, various subsidiary equations and formulas are deferred to appen- 
dices |A| to [|. In the rest of this chapter, we use the units of c = G = 1. 

§2. Teukolsky formalism 

In terms of the conventional Boyer-Lindquist coordinates, the metric of a Kerr 
black hole is expressed as 



A 



sm 



2 , 



ds = — — (dt — asm 9 dip) -\ — — (r + a )d(p — adt 
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+^dr 2 + Ed9 2 , (2-1) 

where S = r 2 + a 2 cos 2 9 and A = r 2 — 2Mr + a 2 . In the Teukolsky formalismE! , 
the gravitational perturbations of a Kerr black hole are described by a Newman- 
Penrose quantity ip4 = —C a p 1 sn a fh /3 n' y fh s , where C a p^s is the Weyl tensor, n a = 
((r 2 + a 2 ),-A,0,a)/(2U) and m a = (ia sin 0,0, l,i/sin6»)/( v / 2(r + ia cos 0)). 
We decompose tp4 into Fourier-harmonic components according to 

(r - mcos0)V4 = E/ due-*** 4 ™' - 2 ^ m (0)/2 /Tnw (r), (2-2) 

im 

The radial function Ri mul and the angular function s S£ m (9) satisfy the Teukolsky 
equations with s = — 2 as 

2 d ( 1 dR hl 



A 2 — [-^—^ ) ~ V(r)R, mw = T imul , (2-3) 
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The potential V(r) is given by 

K 2 + Ai{r - M)K o . 
V(r) = y — J — + 8ioor + \, (2-5) 

where K = (r 2 + a 2 )to — ma and A is the eigenvalue of _2<S^. The angular function 
s Se m (9) is the spin-weighted spheroidal harmonic which may be normalized as 



/ 

Jo 



- 2 Sem\ 2 sm6de= 1. (2-6) 
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The source term T£ maJ is specified later. Here we only mention that for orbits of our 
interest, which are bounded, Ti muJ has support in a compact range of r. 

We define two kinds of homogeneous solutions of the radial Teukolsky equation: 

r Dtrans/\2 -ikr* f nrr ^r, 

^ 1 r^B^e^* + r-^g^e"^* for r +oo, 1 J 



jjup ^IC^+^C^* for r^r + , 
£mw I C%^ s r 3 e^ r * for r +oo 

where k = oo — ma/2Mr + and r* is the tortoise coordinate defined by 

air* 



(2-8) 



"* = / ——dr 
J dr 



2Mr + , r — r+ 2Mr_ , r — r_ , . 

= rH —In In , (2-9) 

r + - r_ 2M r+ - r_ 2M v y 

where r± = M ± \/M 2 — a 2 , and for definiteness, we have fixed the integration 
constant. 

We solve the radial Teukolsky equation by using the Green function method. A 
solution of the Teukolsky equation which has purely outgoing property at infinity 
and has purely ingoing property at the horizon is given by 

Remcu = T7? \ R £mu> f ^'CJfm^ -2 + Rf miJ j [ dr' R^JT^A' 2 \ , 

VVemw \ Jr+ Jr 



(2-10) 



where the Wronskian Wn mw is given by 

W^ = 2iuCtZBZ w . (2-11) 
Then, the solution has an asymptotic property at the horizon as 

otrans a2„— ikr* 

d„^, Zi e 

~<trans p 

The solution at infinity is also expressed as 



otrans a2 — ikr* poo 

RemUr - r + ) - 2? ^ transBinc / dr'RZ^A- 2 ee Zf mul A 2 e^\ (2-12) 
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Here and in the following sections except for section 12, we focus on the gravita- 
tional waves emitted to infinity. Hence will be simply denoted as Z£ mu} . The 
gravitational waves absorbed into the black hole horizon will be treated separately 



in section 12. 



Now let us discuss the general form of the source term Ti muJ . It is given by 

aau) 
lOc, 



where 



with 



Tfyw = 4 J dQdtp-^-\B' 2 + B '^)e- im ^=^, (2-14) 



B' 2 = -±p 8 pL^[p- 4 L (p- 2 p- l T nn )] 

- ^pSpA'L^lp-yj+ip-Zp-lA-'Trnn)], 

B' 2 * = --/-pA 2 J^[p- A J^p- 2 -pT imi )\ 

' p s pA 2 J + [p-VA- l L^(p- 2 p- 2 T mn )}, (2-15) 



2V2 



p = (r — ia cos I 



m 



L s = de H — : — - — aoosinO + scotO, (2-16) 



smi 



J + = d r + iK/A . 



In the above, T nn , T„ n and T^m are the tetrad components of the energy momentum 
tensor (T nn = T^ u n fi n v etc.), and the bar denotes the complex conjugation. 

We consider T^ u of a monopole particle of mass p. The case of a spinning particle 
will be discussed in section |ll| separately. The energy momentum tensor takes the 
form, 

;/ d dz v 

I - ■ - S (r-r(t))6(9-e(t))S( V - V (t)). (2-17) 



rpfMV 



EsinOdt/dr dr dr 

where = (t,r(t),9(t),(p(tfj is a geodesic trajectory and r = r(t) is the proper 
time along the geodesic. The geodesic equations in Kerr geometry are given by 



de 

dt 

dr 



C — cos 2 Q\ a 2 (l - E 2 ) + 
I 



aE 



aE 



sin 



sin 



A 

a sin 2 



sin u 

2 1 „2' 



+ -[E{r 2 + a 



1/2 
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0(6), 
= <P, 
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r 2 + a 2 



A 



E(r 2 + a 2 ) 



T, 



(2-18) 
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where E, l z and C are the energy, the ^-component of the angular momentum and 
the Carter constant of a test particle, respectivelytll E = r 2 + a 2 cos 2 9 and 



R = [E(r 2 + a 2 ) - al z ] 2 - A[(Ea - l z ) 2 + r 2 + C}. 



(2-19) 



Using Eq. (2-18), the tetrad components of the energy momentum tensor are 
expressed as 



T nn = lx—J(r - r(t))5(9 - 9(t))6(<p - <p(t)), 
sm9 

Trnn = ^^S(r - r(t))8(6 - 0(t))% - <p(t)), 
sin 

Tjnm = A*^J<5(r - r(t))5(9 - 9(t))5( V - <p(t)), 
sin u 



(2-20) 



where 



Cn n 
Cmn 

Cmm 



E {r 2 + a 2 )-al z + E^- 
or 



P 



2y/2S 2 i 

r? 

i sin 9 { aE 



E(r 2 + a 2 ) -al z + E 



dr 



i sin 9 ( aE 



sin 



2Et 



sin 



(2-21) 



and t = dt/dr. Substituting Eq. (2-15) into Eq. fl2-14| ) and performing integration 
by part, we obtain 

T tmu = ^L [°° dt [ d 9e iujt ~ im ^ 

y2lT J— oo J 

-ht{p^Lt(p 3 S)}c nn p- 2 p- l 5(r - r{t))5{9 - 8{t)) 
+ (h\s + ia(p - p) sin OS) J+{C- inP - 2 -p- 2 A-H(r - r(t))5(9 - 9(t))} 



+ ^Lt{p 3 S(p 2 p- 4 ), r }C- in Ap- 2 p- 2 5(r - r(t))6(9 - 9(t)) 



^p 3 A 2 SJ + {p- 4 J + (pp- 2 CTnrnS(r - r(t))5(9 - 9(t)))} 



(2-22) 



where 



L\ = da 



rn 



sm( 



+ auj sin 9 + s cot I 



and S denotes _2<S^(0) for simplicity. 



(2-23) 



*' These constants of motion are those measured in units of /i. That is, if expressed in the 
standard units, E, l z and C in Eq. (2T8) are to be replaced with E/fi, l z /[t and C/fj, 2 , respectively. 
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We further rewrite Eq. (2-22) as 



<-lmu — P 



dte itJt ~ imtp ^ A 2 



(A nn + AfnnQ + A m mo)S(r ~ r ( t )) 



+ {(Ajnn! + A mm i)S(r - r(t))} _ + {A^^bir - r(t))} 



, (2-24) 



where 



A nn Q 
AfnnO 

A m m 

Afnn 1 

4 , 

m 1 

A m m 2 



2vrZ\ 2 
2 



C mn p~ 3 \(L+s) {^ + P - -T> 



,K 



-a sin 05— (p — p) 



P pC mm S 



K 



AJ :T A 2 



K 2 K 
+ 2ip- 



V^A 



p 3 Crnn[LtS + iasin9(p- p)S], 



2 _ 3 _^ (.K 

^ P pCrnmbll-^ + p 



P 3 pC mm S. 
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Inserting Eq. (2-24) into Eq. (2-13), we obtain Zi miAj as 



J lmuj 



o,7 , Dine 



dte iwt-im V (t) Wir 



where 



trim) 



^£rnuj{AnnO + A mn Q + ^4 mm o} 



j pin j2 pin 

aIx tmoj J A , 4 \ , " n <mu a 

'X-tt-m nl T ml/ T , 9 m 2 



dr 



dr 2 



r=r(t) 



(2-25) 



(2-26) 



(2-27) 



In this paper, we focus on orbits which are either circular (with or without 
inclination) or eccentric but confined on the equatorial plane. In either case, the 
frequency spectrum of Te muJ becomes discrete. Accordingly Zt muJ in Eq. ( [2-12 ) or 
( 2-13; ) takes the form, 

Z(,mui 

5(u - Un)Z irmjJ . (2-28) 



(2-29) 



Then, in particular, tp4 at r oo is obtained from Eq. ( [2 -21) as 

Imn v 
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At infinity, ^4 is related to the two independent modes of gravitational waves h-\ 
and h x as 

1 



^4 = T^H 



(2-30) 



From Eqs,( 2~2^ ) and ( 2-30| ), the luminosity averaged over t At, where At is 
the characteristic time scale of the orbital motion (e.g., a period between the two 
consecutive apastrons), is given by 



dE\ 



E 



J tmu) n 



E 

£,m,n 



(2-31) 



Imn 



£,rn,n 

In the same way, the time-averaged angular momentum flux is given by 

2 

~~dt 



'cLh 
t dt 



E ■ 

£,m,n 



E 



dim 



E 



(2-32) 
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§3. Post-Newtonian expansion of the ingoing wave solutions 

We consider the case when a test particle of mass ji is in an orbit which is 
nearly circular around a Kerr black hole of mass M ^> fj, and describe a method to 
calculate the ingoing wave Teukolsky functions R^^ which are necessary to evaluate 
the 4PN formulas for the gravitational waves energy and angular momentum fluxes 
emitted to infinity. In the Schwarzschild case, we shall derive the 5.5PN luminosity 
formula in section |j. A method to calculate the ingoing wave solutions in this case 
is separately discussed in Appendix [D| because it is considerably more complicated 
than the method explained in this section. 

Using non-dimensional variables in the Teukolsky equation, we can see that the 
Teukolsky equation is expressed in terms of three basic variables, z = ujr, e = 2Mw 
and au = qe/2 where q = a/M. In order to calculate the gravitational waves induced 
by a particle, we need to know the explicit form of the source terms T^ maJ (r). They 
will be given in the proceeding sections for specified orbits. Here it is sufficient to note 
that they have support only around r = ro where ro is the orbital radius for a circular 
orbit or the mean radius in the case of an eccentric orbit (with small eccentricity). 
Hence from Eq. ( [2T3 ), what we need to know are the ingoing wave functions R m (r) 



around r = ro, and their incident amplitudes B mc . Note that we do not need the 
transmission amplitudes B tTlins to evaluate the gravitational waves at infinity. This 
fact considerably simplifies the calculations. Since we treat a test particle in a bound 
orbit which is nearly circular, the contribution of u) to the Teukolsky functions comes 
from to ~ mQ^, where Q v ~ (M/rg) 1 / 2 is the orbital angular frequency. We will 
evaluate R m by setting three basic variables to be z = u>ro ~ m(M /r^) 1 / 2 ~ v, 
e ~ 2m(M/ro) 3 / 2 ~ v 3 , and auj ~ gm(M/ro) 3//2 ~ v 3 . Here, we have introduced 
a parameter v = (M/ro) 1 / 2 which represents the magnitude of the orbital velocity. 
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We assume that v is much smaller than the velocity of light; »Cl. Consequently, 
we also assume that e -C v <C 1. This relation is the basic assumption in obtaining 
the homogeneous solutions below. 

Now we calculate the ingoing wave solutions which are necessary to calculate 
the luminosity to 0(v s ) beyoncLthe lowest order for the Kerr case. The method is 
mainly based on Shibata et al. ^ and Tagoshi et _al.ll3). An extension to 0(v 11 ) 
calculations done by Tanaka, Tagoshi and Sasaki© for the Schwarzschild case is 
given in Appendix 

First, we discuss the angular solutions. The angular solutions are the spin- 
weighted spheroidal harmonics. The angular equation (2-4) contains only one small 
parameter au. It is straightforward to calculate the spin- weighted spheroidal har- 
monic -2<SV m and its eigenvalue A by expandingthe solution in power of au. It can 
be done by the usual perturbation method^'lill'Ej). It is also possible to obtain 
them by using an expansion by means of Jacobi functions qJ*. The method and the 
results are given explicitly in Appendix [A|. Here we only show the eigenvalue A which 
is used to calculate the radial functions. The eigenvalue A is given by 



A = A + acuXi + a 2 uj 2 X 2 + 0((aiv) 3 



(3-1) 



where 



Ao 
Ai 



with 



-2m 



+ l)-2 = (£-l)(£ + 2), 
£(£ + 1) +4 



+ D 



-2(1 + Dic^f + 2£(c^) 2 + - - ^ + 4)d-3)(l 2 +l-3^) 



■'Im 



Hm 



t+l) 
2 



+ 3)(£-l)(£ + m + l)(£-m + l)] 1 / 2 



(2£ + l)(2^ + 3) 
(£ + 2)(£ - 2)(£ + m){£ - m) 



(2£ + l)(2£-l) 



1/2 



(3-3) 



Next we consider the homogeneous solution R m . We assume u > below. The 
solution for uj < can be obtained from the one for u > by using the symmetry of 
the homogeneous Teukolsky equation which implies Re iin ,uj = Re,-m,-u>- Here, we do 
not treat the Teukolsky equation directly. Instead, we transform the homogeneous 
Teukolsky equation to the Sasaki-Nakamura equationEif, which is given by 



d 2 



_dr* 2 

The function F(r) is given by 



F{r)±-U{r) 



X, 



imuj 



0. 



F(r) 



A 



r] r 2 + a 2 



(3-4) 



(3-5) 
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where 
with 



7] = c + ci/r + c 2 /r 2 + c 3 /r 3 + c 4 /r 4 , 



c = -YZiwM + A(A + 2) - 12au;(auj - m) 
ci= 8ia[3atu — \{aoj — m)], 
c 2 = -24iaM(auj - m) + 12a 2 [1 - 2(aw 
c 3 = 24m 3 (ooj - m) - 2AM a 2 , 



(3-6) 



c 4 = 12a 4 . 
The function J7(r) is given by 



U(r) 



AUi 



(r 2 + a 2 ) 2 



+ G 2 + 



r 2 + a 2 



(3-7) 
(3-8) 



where 



G 
Ui 
a 



2(r - M) 



r 2 + a 2 



+ 



A 2 r 



(r 2 +a 2 ) 2: 

A. 



2a + ^-) -^(n+^i- 



if/3 , 6Z\ 

-?^2" + MK,r + A + — , 

2/y\ 



/9 = 2/A(-ifiT + r-M- ^J. (3-9) 

When we set a = 0, this transformation becomes the Chandrasekhar transformation 
c3 for the Schwarzschild black hole. The Sasaki-Nakamura equation was originally 
introduced, for the inhomogeneous case, to make the potential term short-ranged and 
to make the source term well-behaved at infinity. It is not necessary to perform this 
transformation in this case, since we are interested only in bound orbits. Nevertheless 
we choose to do this because the lowest order solution becomes the spherical Bessel 
function and we can apply the post-Newtonian expansion techniques developed for 
the Schwarzschild case by PoissonEP and SasakilUP. 
The relation between i?^ mw and Xg muj is 



R, 



£mui 



1 If ±PA P 



V 



A 



A' 



(3-10) 



where ximw 



Xg muJ A/ {r 2 + a 2 ) 1 / 2 . Conversely, we can express Xg mu} in terms of 



X £mu) = (r 2 + a 2 ) 1 / 2 r 2 J_J_ 



I/? 



(3-11) 



where J_ = (d/dr) — i(K/A). Then the asymptotic behavior of the ingoing-wave 
solution X m which corresponds to Eq. fl2-7|) is 



y-in 



U: 



trans p —ikr* 



+ A 



imuj 



for r* 
for r* 



oo 



-oo. 



(3-12) 
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The coefficient ,4 inc , A rei and A tr&DS are respectively related to B mc , B ret and B 



inc orcf 



itrans 



defined in Eq. (|2 • 7| ) , by 



nine 


_ /line 


(3-13) 


rjref 


4^ 2 ref 

CO 


(3-14) 


73 trans 


1 /i trans 


(3-15) 
(3-16) 



where cq is given in Eq. (3-7) and 

d lmuj = y / 2Mr + [(8 - 2AiMoj - VoM 2 u?)r\ 

+(12iam - 16M + 16amMw + 24iM 2 uj)r + - 4a 2 m 2 - YliamM + 8M 2 ]. 

Now we introduce the variable z* as 



z + e 

_z + — Z- 
ur* + elne , 



ln(z — z + ) 



z_ 



Z+ — Z- 



■ ln(z — Z-) 



(3-17) 



where z = ur and z± = ojr±. To solve for X m , we set 



Xf mij0 = V z 2 + a 2 u 2 €e m (z) exp (-icp(z)) 



where 



ar u) 

A 



z — z 



In- 



z — z 4 



Z — Z- 



(3-18) 



(3-19) 



With this choice of the phase function, the ingoing wave boundary condition at 
horizon reduces to that £g m is regular and finite at z = z + . 

Inserting Eq. (3-18) into Eq. fl3-4| ) and expanding it in powers of e = 2Mw, we 
obtain 



£ (0) [Him] = 6^ {1) [&n] + eQ W fem] + <?Q W fern] + e 3 Q ^ [6m] + [6m] + 0(e b 

(3-20) 

where L^°\ and are differential operators given by 



c2n( 2 )f 



,3n(3)[ 



, 4 n( 4 )[ 



L (0) = ^_ + 2_d 

<iz 2 z cfe V z 2 
L(1) = 1^ + / 1 + 2 ! W _ 

z (iz 2 V z 2 z J dz 



4 t 1 

z 3 z 2 z 



Q 



(1) 



iq\\ d Aimq 



2z 2 dz 1(1 + l)z 3 ' 



(3-21) 
(3-22) 
(3-23) 
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The formulas for , and are very complicated, and they are given ex- 
plicitly in Appendix |B|. Note that, when we set a = 0, all vanish. 
By expanding ^£ m in terms of e as 

oo 

^m = J2e n ^(z), (3-24) 

n=0 

we obtain from Eq. ( |3-20| ) the iterative equations, 

L^$>]=z-*WM, (3-25) 

where 

= 0, (3-26) 
W&=*(LW[&]+QW[$h), (3-27) 

= z 2 K£] + Q W K&] + Q (2) [&]) , (3-28) 
Wfi = z 2 m + Q W i&] + Q (2) K&] + Q (3) i&) , (3-29) 
= - 2 (* (1) K2] + Q (1) KS] + Q^lSl + Q (3) [a + Q (4) [^])(-3-30) 

As e = 2GMu if we recover G, the above expansion corresponds to the post- 
Minkowski expansion of the vacuum Einstein equations. 



The iterative equations (|3-25|) have been obtained by expanding Eq. (|3-4j) in 
powers of e by regarding z = cor as the independent variable. Since the horizon is 
at z = z+ = 0(e), this procedure implicitly assumes that e <C z. Consequently, we 
cannot apply the above expansion near the horizon where the ingoing wave boundary 
condition is to be imposed. To implement the boundary condition correctly, we have 
to consider a series solution of which is valid near the horizon as well as in the 
range e <C z and match it to the series solution of the form (3-24). Recenth^this 



matching problem has been rigorously solved by Mano, Suzuki, and Takasugic^) for 
the original Teukolsky equation. However, for our present purpose, it is sufficient 
to resort to a simple power-counting argument, by which it is possible to implement 

(n) 

the boundary condition of ^ m at the horizon to the behaviors of Q m at z > e for 
n <2i (for n < 21 + 1 in the Schwarzschild case; see Appendix [D]). 

Since the ingoing wave boundary condition is that £g m is regular at horizon, if 
we introduce an independent variable i := (z-z+)/ewe can expand £g m near the 
horizon as 

oo 

&m = £ e %t } (x>> (3-31) 
n=0 

This means we have 

oo 

6 m (0)=^ } (0)5>"C„, (3-32) 

n=0 

where C n = ^(0)/Q^ (0). In other words, & m (0) should have a well-defined limit 

for e — > except for the overall normalization factor C/^?(0). Keeping this property 
in mind, let us consider the boundary conditions for Eqs. (3-26) - (3-30). 



14 



The general solution to Eq. (3-26) is immediately obtained as 

where jg and ng are the spherical Bessel functions. The coefficients and (5^ 
are to be determined by the boundary condition. For convenience, we normalize the 
solution X m such that the incident amplitude A mc is of order unity. Then both ot^ 

and (3g^ must be of order unity. Since jg(z) ~ z l and ng{z) ~ z - ^" 1 for z<l, the 
latter is 0(e 2 ^ +1 ) larger than the former near the horizon where z = O(e). Hence 
Eq. ( 3-32j ) implies we should set = 0. As for the value of a^, since it only 



contributes to the overall normalization of X m , we set al^ = 1 for convenience. 

Inspection of Eqs. (3-27) - (3-30) reveals that the solution behaves as z e ~~ n 
plus the homogeneous solution ocj^jg + fig^ng for z — ► 0. As for (n > 1), they 
simply contribute to renormalizations of otgl,- Hence we put them to zero. As for 

in) 

Pl m , from the same argument as given above, we find they may become non-zero 
only for n > 2£ + 1. Since I > 2 and e = 0(v 3 ), this implies that the near zone 
contribution of ng ~ z~^~ l , which is 0(t>~( 2 ^ +1 )) greater than the lowest order term 
jg, to the gravitational waves emitted to infinity may arise only at 0(v 10 ) beyond 
the quadrupole order. Since the post-Newtonian corrections we shall consider for the 
Kerr case are those up to 0(v 8 ), we set = and solve the iterative equations 
( 3-25| ) to 0(e 4 ) with the boundary conditions that ~ z e ~ n at z — > 0. We note 



that, in the Schwarzschild case which is discussed separately in Appendix [S|, these 
boundary conditions turn out to be appropriate for n < 2£ + 1; i.e., up to one power 
of e higher than the Kerr case. 

To calculate f° r w > 1, we rewrite Eqs. (3-27) - (3-30) in the indefinite 
integral form by using the spherical Bessel functions as 



= /" dzjeWiS ~ k f dzngW^ . (3-34) 



The calculation is straightforward but tedious. All the formulas which are needed 

(n) 

to calculate the above integration to obtain Q m for n < 2 are shown in Appendix of 
Sasaki 0) . They are recapitulated in an alternative way in Appendix |D|. Using those 
formulas, we have for n = 1, 

(1 ) _ (l-l)(l + 3) . 



2(£ + l)(2£ + l) J ^ x \2£{2£ + l) £{£-! 




+Re,ojo + Y] ( — H ^— r ) Re,mjm - 2D n l i + ijg In z 

m=l 

img / £2 + 4 \ . imq ( (l+l) 2 + 4 \ . 



2 \ y £ 2 (2£ + l) y / J 2 ^(£+ 1) 2 (2£ + 1) 

where D? J is an extension of the spherical Bessel function defined in subsection D.2.1 
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of Appendix O; 



D? = - \j t Si(2z) - n t (Ci(2*) - 7 - In 2z)\ , (3-36) 



where 7 = 0.5772- •• is the Euler constant, Ci(x) = — f£° dtcost/t and Si(x) = 
Jq dtsint/t, and R m ,k is a polynomial of the inverse power of z defined by 

Rm,k = % {rimjk jmlT'k) 

[(m-fc-l)/2] ( m _ /. _ I _ r \\ f ( m+ 1 ^ /9 vm-fc-l-2r 

= - E (-1)- — K-^ — M-) (3-3?) 

t^o r! (m-k — 1 — 2r)! rf/c + f + rj W 

for m > k and it is defined by 

Rm,k = —Rk,m- (3-38) 

for m < k. 

Next we consider Q 2 ^. From Eqs.( |3-34| ) and (3-35), we obtain as 

Z?2 = fP + 4 2) + 4%) i (3-39) 
where fj, 2 ' and gP' are the real and imaginary parts of in the Schwarzschild 

(2) 

limit, respectively, and fc)~(<z) is the correction term due to non- vanishing q = a/M. 
For £ = 2 and ^ = 3, /j 2 ^ are given by 



(2) 389 . 113 . 1 . j 

2 = -70^ J0 - 420l Jl + 7; J3 + 4D2 
5 „n,o 10 6 „„,■ 107 



3z 2 3 1 z 105 ~ 3 
107 . 1 , 2 

-210 J2lnZ "2 J2(ln2) ' 



,( 2 ) 1 . 323 . 5065 . / 1031 445 \ . 65 65 
" ~ ' 1 + Wz J2 ~ 294? J1 " l588l + T4?J J ° + 6?"° " 6? 



3 n ni , 13 n j 9 nj 30 nj - 13 nj 
z 3 + 3 2 + z 1 + ^ V ° "21 - 4 
+4 J D3^-^ j3 lnz-- J 3(lnz) 2 , (3-40) 

where -D™ n-? is defined in Appendix^, Eq. (D-20). As explained in subsection D.l 
of Appendix [D], the term g^p is given for any £ as 

gf ] = - l -j f , + fP\nz. (3-41) 

(2) 

The term k^(q) is given for £ = 2 and ^ = 3 by 

( 2 ) 191 i . m 2 q 2 j mqji 
k\' = mq in 

2m 18Q *JU 3Q 1Q 



16 



68 i 
~63~ 
11 i 



73 m 2 » 2 



mgj 2 



392 



J2 



( 3i 
+mq — — - 



q 2 h 

392 
13 h 
90 



1764 
71 m 2 q 2 j\ 
8820 

In z + imq 



<T . . 7mqj 3 
32 + 



+ 



180 
13 i 
~6 



% o X 9 9 

72<7 is + ^m g j3 



mqn\ 



4d? 



13 n j 
45 3 



(3-42) 



,(2) 3527 i . 2m 2 q 2 ji mqj2 hi 2 . hi 2 2 . 

Kn'= mqii g 72 H m q 70 

3m 84Q hji 315 36 5Q4 «i ^ 226g w .« 

379i . q 2 j 3 7m 2 g 2 j 3 3mqj A i 2 . i 2 2 . 

m g 73 g 74 H m g 74 

360 360 720 160 140 1120 

97 1 g 2 j 5 17 m 2 q 2 j 5 103 1 25 1 

+ 504fl m9j5 - 360 _— 5040 4 ^mgn + — mgn 2 

- f 13mgj2 + In * + img f Zl? - 1 D ?A . (343) 

V 126 56 J y V 63 2 28 4 J v 7 

As noted previously, the source term T^ muJ has support only around r = r , 
hence around z = euro = 0{v). Therefore, to evaluate the source integral, we only 
need X m at z = 0(v) <C 1, apart from the value of the incident amplitude A mc . 
Hence the post-Newtonian expansion of X m corresponds to the expansion not only 
in terms of e = 0(v 3 ), but also of z by assuming e <C z <C 1. In order to evaluate the 
gravitational wave luminosity to 0(v 8 ) beyond the leading order, we must calculate 
the series expansion of q"^ in powers of z for n < 6 — I for each 2 < I < 6. This 
follows from a simple power counting. The leading order contribution of the l-th. pole 
is 0(v 2 ^~ 2 ^) smaller than that of the quadrupole, while the n-th post-Minkowski 
terms are 0(e n z~ n ) = 0(v 2n ) relative to the lowest order terms in the near-zone. 
Hence the leading term of £^ contributes at 0(v 2 ^~ 2 ^ +2n ) and 0(v 8 ) is attained 
for I + n = 6 (see Appendix C of Ref. [jX 



(n) 

To evaluate A mc , we need to know the asymptotic behavior of <Q m at infinity. 
Since the accuracy of A mc we need is 0(e 2 ), we do not have to calculate and 
in closed analytic form. We need only the series expansion formulas for and 



£(m aroun< i z = 0, which are easily obtained from Eq. (3-34). This is also true for 
£4^ for I = 4. Inserting £^ into Eq. (|348 ) and expanding it by z and e assuming 
e < 2 < 1, we obtain 

( 3 ) -g 2 i 3 -i Imq i 2 2 mg 3 m 3 g 3 

£9™ = m g H m q H 

S2m 30z 30z y 30 180 60 * 36 90 

m q In 2 i 2 2 

m g In 2 

30 

100637 i q 2 17 m 2 g 2 83 i , 

H m g 1 1 m g 

441000 180 1134 5880 

, , In z 106 i . i , 2 

■ m q H mq mz mq (mz) 

13230 15 1575 30 v ; 
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f (4) 

S2m 
e(3) 

t(2) 
?4m 



„4 



80 z 
—i 
1260 
/ 1 



mq + 



1 2 

m g 



V1764 15120 



1890 
11 i 

mq + 



1260 



m g 



- 3780 m3g3 + O( ^ 
19 m 2 g 2 ' 



10584 105840 



z z + 0(z A 



(3-44) 
(3-45) 

(3-46) 

(3-47) 



Inserting these £^ into Eq. ( |3-18|) and expanding the result in terms of e = 2Mu>, 
we obtain the post-Newtonian expansion of X m . The transformation from X m to 
R m is done by using Eq. ( 3-10 ). 



Next, we consider A mc to O(e ). Using the relation ~ —ji-i 



\l+n 



n2n-e 



at z ~ oo, etc., we obtain the asymptotic behaviors of an d dm at z ~ cxd as 



(2) 
£m 



+ ln 2 ~ ( ln 2 ) 2 ) ^ + («£m - rim 111 Z ) H t 



+ipflje^z + i(q£l 



where 



D (i) 

y^m 



7T 

~2' 
1 

2 

£-1 



In z)ri£ ln z 



l)(£ + 3) 



+ 1 



ln2 



2imq 
¥{£ + !) 



fcr fc 



7, 



(3-48) 

(3-49) 

(3-50) 
(3-51) 

(3-52) 



for any I and 



J2) 
F2m 



(2) 
92m 

(2) 
P3m 



fl (2) 
°3m 



4577 _ 
210 ~ 

-7 ln2 



7^ 5tt 2 

— H 

2 24 



— 7 m g + 
18 ' 



457 In 2 
210 



i (ln2) 2 

— mo ln 2 . 

18 2 



-457 7T 7 7r 5mg j 

1 1 1 rrnrq 

420 2 36 36 y 



52 7 
~2T 



7 2 5^ 2 

— H 

2 24 



72 



-7 ln 2 mo ln 2 - 

' 72 y 

-267r 77T 67mg 



jmq + 
(ln2) 2 



72 
52 ln2 
21 



i 22 vr ln2 

H m o H 

324 y 2 



(3-53) 
(3-54) 



21 



+ 



H mir q ^ o 

1440 144 H 360 



12960 



2 2 i 

m g + 



(3-55) 

7T In 2 . 
-^6) 



Then noting that exp(— i 
is expressed as 



j exp(— — z)) at z ~ oo, the asymptotic form of X 1 - 



18 



(i)> 



T (2)> 



(3-57) 



where fo^ and np are the spherical Hankel functions of the first and second kinds, 
respectively, which are given by 

hf> = k + in, -> , hf = u ~ ini - (-l) m£ ^- 

From these equations, noting oor* = z* — elne, we obtain 



(3-58) 



^inc _ _^+l g — ieln< 



l + e(p{2 + ^)+e 2 (^ + ^) + - 



2^(2) 



,(2)n 



(3-59) 



The corresponding incident amplitude B mc for the Teukolsky function is obtained 
from Eq. (3-13). 

§4. Gravitational waves to 0(v n ) in Schwarzschild case 



In this section we consider a circular orbit around a Schwarzschild black hole 
and derive the 5.5PN formula for the energy flux emitted to infinity. In this case, we 
can take the orbit to lie on the equatorial plane {6 = vr/2) without loss of generality. 
Then E and l z are given by setting R(ro) = dR/dr(ro) = where R is given by 
Eq. (HID . This gives 



E = (r - 2M)/ v /r (r -3M), l z = y/M^/fi - 3M/r , 



(4-1) 



where ro is the orbital radius. The angular frequency is given by Q v = {M/tq) 
Defining s b lm by 



3\l/2 



7T 



Ohm = g - + + 2 )] 1/2 O^m ( -• 

.x6, m = p - 1)(£ + 2)] 1 / 2 „ x Y lm (~, o) ± 

\2 J r 

-•film = -2^m [ ~2>®) Iz^Vi 



Er 



r - 2M 



(4-2) 



where s Yf m (9, <p) are the spin-weighted spherical harmonics 0, Z£ muJ is found to take 
the form 

Zimcu = Z im 5(uj - Tnf2tp), (4-3) 

where 



7T 



2i^b em (l + Jwrg/(r - 2M)) 



+i_2&< m wro(l - 2M/r )- 2 (\ - M/r + ^r Q J 
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+ 



i-ibe m --2 h m (l + iujr^/(ro - 2M)) 



roRZ (ro) 



+^-2k m r 2 Q RZ (tq) 



(4-4) 



where prime denotes the derivative with respect to the radial coordinate r. In terms 
of the amplitudes Zi m , the gravitational wave luminosity at infinity is given by 



ira I 



(4-5) 



where uj = mfl^. Since the dominant frequency of the gravitational waves at infinity 
is 2f2,n, an observationally relevant post-Newtonian parameter is x = (Mi?^) 1 / 3 . We 
mention that our post-Newtonian expansion parameter is defined by v := (M/ro) 1//2 . 
In the case of a circular orbit around a Schwarzschild black hole, however, we have 
v = x. Hence the parameter v is directly related to the observable frequency in the 
present case. 

Following the method given in section pi instead of directly calculating R 1 ^ from 
the homogeneous Teukolsky equation, we calculate the corresponding Regge-Wheeler 
function X"^ first and then transform it to R"^- The homogeneous Regge-Wheeler 
equation, which is given by setting a = in Eq. (|3-4j), takes the formLiP, 



d 2 
dr* 2 



V(r) 



XUr) 



0, 



where 



V(r) 



r J V 

The transformation ( |3 • 1 0| ) reduces toi 2 



2M \ ( £(£ + 1) _ 6M \ 



(4-6) 



(4-7) 



Aid \r 2 (d . ... 

R£ ™ = ^{d^ + w )A{-d^ + w]rXe " 



where cq, defined in Eq. (3-7), reduces to cq 
inverse transformation ([3-1 1|) reduces to 



(4-8) 

!)£(£ + !)(£ + 2)- YliMuj. The 



X, 



r 5 / d 
~A \dr* 



iuj 



d 



A \dr* 



IUJ 



R 



(4-9) 



The asymptotic forms of X m is the same as given in Eq. ( 3-12 ) except that now we 
have k = oj. The coefficients A mc , A rei and A trBjUS are also respectively related to 
B' mc , B TC{ and 5 trans as before. (See Eqs. (3-13), (3-14) and (3-15).) Note that the 
coefficient that appears in Eq. (3-15) now reduces to 

dimoj = 16(1 - 2iMw)(l - 4*Afw). (4-10) 

Corresponding to Eq. (3-17), we introduce the variable z* = z + eln(z — e). Then 



Eq. (3-18j) reduces to 



-i(z*—z) 



-it ln(2— e) 



(4-11) 
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and Eq. ( 3-20| ) becomes L^ )^] = eL' 1 ' [&] . Thus Eq. ( 3-25| ) simplifies considerably 
to become 

L(°)[^ n) ] =L«[cJ n_1) ]. (4-12) 
It may be worthwhile to note that the left hand side can be expressed concisely as 



dz 



Id/ 

z 3 dz V 



e iZ z 2 &~ 1] 



(4-13) 



The calculations to 0(e 2 ) are already done in section |3|. When we consider the 
gravitational wave luminosity to O^v 11 ), we need to calculate ^4 mc to 0(e 3 ) for £ = 2 



and 3 and to 0(e) and for £ = 4. Thus we need the closed analytic forms of £j for 

(2) 

£ = 2 and 3 and Q . The latter can be obtained in the same way as in the previous 
section. The procedure to obtain ^ is explained in detail in Appendix 



The real parts of £ 



(3) 



f?\ for 



2 and 3 are given as 



( 3)_214F 2i0 [z(lnz)jo] 107 ^ 4 457 D 



105 

107 (Inz) Dl 1 



105 

_ 8D nnnj _ 



630 
+ (In z) 2 D^ 



70 

2D n J 

49 



) "' 2629 Dq J 16949 D% 3 



630 4410 
12^-18 IT" 



2D™ J 



197 j_ 3 2539 j-i 107(lnz)j_i 3(lnz)'j_i 21 ^ 



126 



3780 



349 (In z) ji 9(lnz) 2 ji _ 457 j 3 
140 " 4 1050 



+ 



70 

29 (lnz)j 3 
252 



Qnz) 



12 



100 



(3) = 26F 3 ,o[z(logz)j ] 13 14075 £g| 

/3 21 98 882 

1424 Dl\ 2511 ^ 269 D^ 3 13 (log 2) Iff J 
63 70 + 70 21 

+(log zf - 3£ + 60 + 34 - + 

lo zi ( 

_ 8jD nn n j 75j_4 _ 19j_ 2 _ 65 (logz)j- 2 _ 15(logz) 2 j. 2 _ 2789 j 

3 28 56 14 2 3780 

221 (log z) j _ 17 (log zf j _ 7495 j 2 4867 (log z ) j 2 

84 4 5292 1764 

355(logz) 2 j 2 10963 j 4 15 (log 2:) j 4 3(logz) 2 j 4 4j 6 

84 32340 196 28 1485 ' 1 ' j 



where the definitions of the functions £ff nj etc. are given in Eqs. (D-20) and ( D-24 ) 
of Appendix |5[ The imaginary parts are expressed in terms of and fj, 2 ^ as 

(2) (2) 

given in Eq. (D-9). As for the real part of Q , fl , it is calculated to be 

(2) _ 56 . / 5036 30334 \ . / 35252 30334 14401 \ . 

iA ~ 16o7 J5 + v~33^ + 1155^2 J J4+ v 33^~ 165^ + 3465^ J J3 
/5036 45461 36287 \ /140 5 \ 49 

-\^ + m^ + mm-z) 3l + \—-Wz) n °-^ n2 
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__ r> nj 4- — n nj - — n" 3 4- — n n3 + — - — n nj 

5z 4 30 3 a* 2 - 2 1 - 3 ~ 



+4D nnj 



149 

liO" 
1571 

6930 



10 



1 



105 



210 

T 3 



20 

2; 



1571 



3465 



-D 



j 4 In z - - j 4 (In z 



(4-16) 



Using the analysis given in subsection D.4.2 of Appendix O, the above results 

(3) (2) 

readily give us the asymptotic forms of Q (£ = 2, 3) and Q at z — > 00, from which 
the amplitudes A mc to the required order are calculated. The results are 



^3 



4 mc 

^4 



1 



• e -i e (ln 2e+7 ) exp 



2 107 o / 29 

e 2 7r + e 3 

420 V648 



7T 2 

e 2 +£ 



25 5 2 107, , n , 

1 n H (7 + In 2) 

18 24 210 w ' 



25 
36' 



-7T 



107, . 7T 3 

(7 + ln2)7r 

420 w ; 16 



+ 



e -ie(ln2e+ 7 ) exp 



13 



29 
810 



7T 2/ 169 5 2 13, , A 

1 - e- + e 2 + — 7T 2 + —(7 + n2) 

2 V 72 24 42 w ') 



, / 169 13 , , „ s 7T 3 \ 

+e 7r (7 + ln2)7r + 

I 144 84 w ; 16 ) 



i ie -ie(ln2e+ 7 ) 

2 

x 



cxp 



149 
~60~ 



107 2 C(3)\ 

7T H 

1260 3 J 



+ 



252 3 J 



2 1571 

7T + 



7T { 2 / 22201 
V 7200 



1 - e- + e 
2 



13860 
5 2 1571 

H 7T^ H 

24 6930 



(7 + In 2) +■ 



(4-17) 



The corresponding amplitudes B mc are readily obtained from Eq. (3-13). 

As in the previous section, from Eqs. ( 4T1| ) and ( f4-8|) , it is also straightforward 
to obtain the near-zone post-Newtonian expansion of X m and hence of R m , assuming 
z<l. As discussed there, we need the series expansion formulas for R m for 2(n + 
I - 2) < 11, hence for n < 7 - I for each 2 < I < 7. The resulting R in for 2 < £ < 7 
which are necessary to calculate the luminosity to 0(-u n ) are given in Appendix ^. 

Finally, from Eq. ( f4-5|) , we obtain the luminosity to 0(v 11 ) as 
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where (dE/dt)^ is the Newtonian quadrupole luminosity given by 

'dE\ _ 32/x 2 M 3 _ 32 / fi 
~)n ~ 5rg ~~5\M 

To compare the above result with those obtained previously by the standard post- 
Newtonian method, we note that v = x = (Mi?^) 1 / 3 in the present case. Then 
we ^dpur ij£§ult agrees with the standard post-Newtonian results up to 0(x 5 ) 
g3(),|34|),|35|),gg),g7p,|38|) j n \[ m [^ ^jM <C 1. The contributions to the luminosity from 

individual I modes are given in Appendix [E|. 



§5. Convergence of the post-Newtonian expansion 



Using the results obtained in the previous section, we compare the formula for 
the gravitational wave flux with the corresponding numerical results and investigate 
the accuracy of the post-Newtonian expansion. 

A high precision numerical calculation of gravitational waves from a particle in a 
circular orbit around a Schwarzschild black hole has been performed by Tagoshi and 
Nakamura.H2P Since no assumption was made about the velocity of the test particle, 
their results are correct relativistically in the limit [i <^M. In that work, dE/dt was 
calculated only for £ = 2 ~ 6. Here, for the orbital radius ro < 100M, we calculate 
dE/dt again for all modes of t = 2 ~ 6 and for 1 = 7 with odd m. The estimated 
accuracy of the calculation is about 10 -11 , which turns out to be accurate enough 
for the present purpose 1 _A.s for the radius ro > 100M, we use the data calculated by 
Tagoshi and Nakamurali3) which contain modes from £ = 2 to 6. 

In Figs. 1 and 2, we show the error in the post-Newtonian formulas as a function 
of the orbital radius r. The error of the post-Newtonian formula is defined as 



error = 




where (dE/dt) n and (dE/dt) denote the (n/2)PN formula and the numerical result, 
respectively. In the plot of Fig. 2, only the contributions from I = 2 to 6 are included 
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error 




Fig. 1. The error of the post-Newtonian formulas as functions of the Schwarzschild radial coordinate 
r for 6 < r/M < 100. Contributions from i = 2 to 7 modes are included. 



in both the post-Newtonian formulas and the numerical data. We can see that, at 
small radius less than r ~ 10M, the error of the 1PN and 2.5PN formulas are larger 
than the other formulas. On the other hand, the Newtonian and the 2PN formulas 
are very accurate within this radius. This is because those formulas coincide with 
the exact one accidentally at a radius between 6M and 10M. The error of each 
post-Newtonian formula at the inner most stable circular orbit, r = 6M, becomes 
as follows; 12% (Newtonian), 66% (1PN), 8.6% (1.5PN), 3.4% (2PN), 42% (2.5PN), 
11% (3PN), 5.4% (3.5PN), 17% (4PN), 8.4% (4.5PN), 6.5% (5PN), 4.1% (5.5PN). 
As is expected, the errors of the post-Newtonian formulas decrease almost linearly 
up to r ~ 5000M in a log- log plot. This fact also suggests that the numerical data 
have accuracy of at least ~ 10~ 18 at r ~ 5000M. 

In order to examine exactly to what order the post-Newtonian formulas are 
needed to do accurate estimation of the parameters of a binary, using data from 
laser interferometers, we must evaluate the systematic error produced by incorrect 
templates. However, here we simply calculate the total cycle of gravitational waves 
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error 




100 1000 



Fig. 2. The same figure for 100 < r/M < 5000. Contributions only from I — 2 to 6 are included in 
both the post-Newtonian formulas and in the numerical data. 



from a coalescing binary in a laser interferometer band and evaluate the error pro- 
duced by the post-Newtonian formulas. It has been suggested that whether the error 
in the total cycle is less than unity or not gives a useful guideline to examine the 
accuracy of the post-Newtonian formulas as templates^ (see also Ref. flO|)). 

We ignore the finite mass effect in the post-Newtonian formulas and interpret 
M as the total mass and fx as the reduced mass of the system. The total cycle N 
of gravitational waves from an inspiraling binary is calculated by using the post- 
Newtonian energy loss formula, (dE/dt) n , and the orbital energy formula (dE/dv) n 
which is truncated at n/2PN order as 

^ = rjkm»± (5 .2, 

J Vf 7T \{dE/dt) n \ 

where Vi = (M/ri) 1 / 2 , Vf = (M/rf) 1 / 2 , and and rj are the initial and final 
orbital separations of the binary. We define the relative difference of cycle AN^ 
as AN^ = |iV (n) - i\r( n-1 )|. We adopt r f = 6M and n is the one at which the 



Chapter 1 Black Hole Perturbation 



25 



frequency of wave is 10Hz and which is given by n/M ~ 347(M /M) 2 / 3 . 

The results for typical binary systems are given in Table I. We only show the 

Table I. The relative difference of cycle AN^ for typical coalescing compact binaries. The last 
line shows the cycle calculated by Newtonian quadrupole formula. 



n 


(1.4M ,1.4M©) 


(1OM ,1OM ) 


(1.4M ,1OM ) 


(1.4M ,7OM ) 


2 


356 


54 


216 


212 


3 


228 


60 


208 


296 


4 


11 


5 


15 


31 


5 


12 


7 


20 


53 


6 


11 


8 


22 


75 


7 


1.2 


1.0 


2.6 


10 


8 


0.12 


0.14 


0.3 


2.2 


9 


0.82 


0.80 


1.9 


8.9 


10 


0.09 


0.08 


0.20 


0.87 


11 


0.03 


0.03 


0.07 


0.40 


N w 


16000 


600 


3578 


898 



results_for q = 0. The cases for q ^ are investigated in Shibata et al.ll^P and Tagoshi 
et al.li^P . This table suggests that we need the 3PN~4PN formula to obtain accurate 
wave forms for typical binaries whose total mass are less than 20M & . Although the 
required post-Newtonian order is very high and it has not been achieved yet in 
the standard post-Newtonian analysis, this results show that the post-Newtonian 
approximation is applicable to the inspiral phase of coalescing compact binaries. In 
this sense, we can be optimistic. 

On the other hand, the convergence for the case of neutron star — black hole 
binaries, whose mass is above several ten Mq, is very slow. This is because ti/M 
become smaller for a larger mass black hole, and the higher relativistic correction 
becomes more important. From Table I, we mig ht think that iV (ri) converges at 
n = 11 even for (mi, 7712) = (1.4M0,7OMq). However this is not true. Note that 
Table I shows only the relative difference between the post-Newtonian approximated 
cycles. If we calculate the difference between the post-Newtonian formula and the 
fully relativistic one, we find that the 5.5PN formula is not accurate enough forthe 
case (7711,7712) = (1.4M ,7OM ), as pointed out by Tanaka, Tagoshi and Sasakill3. 

Finally we comment on the initial frequency. The above results are obtained 
by setting the initial frequency to 10Hz. However, it may be difficult to observe 
gravitational wave at this frequency because of the seismic noise. If we set the 
initial frequency higher than 10Hz, the error AN becomes slightly smaller. But 
since this dependence of AN on the initial frequency is very weak, the above results 
are insensitive to the variation of the initial frequency. 

§6. Circular orbit on the equatorial plane around a rotating black hole 



In this section, we consider a circular orbit on the equatorial plane of a Kerr 
black hole and calculate the 4PN luminosity formula. 

We define the orbital radius as r = rn. As in section |4], we have C = 0, and E 
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and l z are determined by R(vq) = and dR/dr 



r=ro " 



as 



E 



1 - 2v 2 + qv 3 



rov(l — 2qv 3 + q 2 v ^ 



(1 - 3v 2 + 2qv 3 y/ 2 ' 2 (1 - 3v 2 + 2qv 3 y/ 2 
where v = (M/ro) 1 / 2 . From these, we can easily obtain (p(t) as 



<p(t) = Q v t ; i? ( 



When a = 0, this becomes Q a 



_ M 1 / 2 r 

^ ~~ 3/2 



1 - qv 3 + gV + 0(v 9 ) 



(6-1) 



(6-2) 



(M/r, 



3\l/2 



The rest of the calculation is almost the same as in section ||. The amplitude of 
the Teukolsky function at infinity is expressed as 
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(6-3) 



where A nn o, etc. are given by Eq. (2-25). 

The total luminosity up to 0(v 8 ) is expressed as 
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(6-4) 



where (dE/dt)^ is the Newtonian quadrupole luminosity, Eq. ( [4- IE ), and the q- 
independent terms are identical to those in Eq. (4-18). 

From an observational point of view, it is more convenient to express the luminos- 
ity in terms of the variable x = (M J?^) 1 / 3 . Using the relation between v = (M/ro) 1 / 2 
and x given by Eq. (6-2), the luminosity is expresses as 
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/ 3597rg | 22667 q 2 | 17q r 



14 



4536 



16 



(6-5) 



where 



(6-6) 



dt / N 5 \M J 

and the q- independent terms are again identical to those in Eq. (4-18) with the 
replacement v — > x. The partial luminosities for individual modes are given in 
Appendix [|]. Tlmspin dependent term at 0(v 3 ) agrees with the standard post- 
Newtonian resultclP. 

Here it is interesting to investigate the origin of some of the spin-dependent 
terms. As an example, we consider the mode i = \m\ = 2. The luminosity from the 
I = \m\ = 2 modes is given by 
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(6-7) 



We can derive some of the spin-dependent terms in the above formula from the 
quadrupole formulae^; dE/dt = (32/5)^ 2 f 4 f?^, where r is the orbital radius of a 
test particle in harmonic coordinates. If multipole moments of the black hole exist, 
the orbital radius changes due to the influence of those moments. The mass and 
mass current multipole moments of a Kerr black hole is given by M; +iSi = M(ia) 1 . 
We can express the orbital frequency of the test particle in harmonic coordinates. 
We find that the dominant effect of the multipole moments of a Kerr black hole to 
dE/dt can be expressed as 
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The terms in this expression agree with the corresponding terms of our result such 
as (— 8/3)qx 3 , 2q 2 x A , (—8/3)q 3 x 7 and q 4 x 8 . Thus, we may interpret the term 2q 2 x 4 
as the effect of the quadrupole moment. The terms (— 8/3)q 3 x 7 and q x 8 are not due 
to a single multipole moment, but to combined effects of the multipole moments. 
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§7. Slightly eccentric orbit around a Schwarzschild black hole 

In this section, we present post-Newtonian formulas of gravitational waves from 
a particle in slightly eccentric orbits around a Schwarzschild black hole. We derive 
the 4PN formulas of the energy and angular momentum fluxes to 0(e 2 ) where e is 
the eccentricity of the orbit. 

The solution of the geodesic equations for slightly eccentric orbits has been given 
by Apostolatos et al.c 3 !. Here we briefly sketch the derivation of it. Since we may 
consider the orbit to lie on the equatorial plane without loss of generality, we put 
6 = ir/2 and C = in the geodesic equations (2-18). Then we define a slightly 
eccentric orbit as follows: First of all, we assume that E and l z are set to be such 
that R(r)/r 4 , which plays the role of an effective potential for the radial motion, 
has the minimum at r = tq and that the maximum value of the orbital radius is at 
r = ro(l + e), where eCl. Thus, the following conditions hold; 

d{R/r A ) , = . = Q R r r= t 1 + e \\ = (7.!) 

or 

From these equations, E and l z are expressed in terms of ro and e as 

2 (l-2v 2 ) 2 v 2 {l-6v 2 ) 2 2v 2 {l-7v 2 ) o 4n 
E 2 = \ 9 J + 5-5 V " T 515 ^ + ^(e 4 ), (7-2) 
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where u = ^/M/tq. For convenience, we also define i7 = u 3 /M_Then expanding the 
geodesic equations in powers of e, the solution is found to beE^ 

r(t) = ro[l + ecos /? r i 

+e 2 {gi(t;)(l - cos i? r i) + gs(u)(l - cos2J2 r i)}] + 0(e 3 ), (7-4) 
ip(t) = Q v t — epi(v) sin Q r t 

+e 2 {P2(v) sin Q r t + v${v) sin 2^ r i} + 0(e 3 ), (7-5) 



where 



Q r = 12(1 - 6?; 2 ) 1/2 , (7-6) 

,/ n 2l n 3(1 - 3w 2 )(l - 8t; 2 ) 

^ = fl[l-/(^]; /(.) - J—JL—J , ( 7 . 7 ) 

. . 1 - 7u 2 , . 1 - llu 2 + 26v 4 , . 

^ ) = T^2' ^ ) = 2(l-6, 2 )(l-2, 2 )' (7 ' 8) 
2(1 -3w 2 ) , , 2(1 - 3f 2 )(l - 7v 2 ) 



(l-2u 2 )vT=fo?' (l-2t> 2 )(l-6t; 2 ) 3 / 2 ' 
_ 5 - 64t> 2 + 250?; 4 - 300?; 6 
PsW ~ 4(1-2?; 2 ) 2 (1-6 V 2 ) 3 / 2 ' ( } 

As is well known, since Q r ^ JX>, the orbit does not close. 
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Now we evaluate the source term of the Teukolsky equation. In the present case, 
A nn o etc. given in Eqs. (2-25) reduce to 
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s£! = L\4Se m (0 = vr/2), 



5S=4^ m (^ = vr/2), 
sQ = S em (9 = ir/2). 

Then noting that the orbits of our interest have the properties, 
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where At r = 2ir/f2 r , Eq. ( ggg ) can be rewritten as 
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u; n = ni? r + m/?^ (n = 0, ±1, ±2, . . .). 



(7-10) 
(7-11) 
(7-12) 
(7-13) 
(7-14) 
(7-15) 



(7-16) 
(7-17) 
(7-18) 
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We see that Zi muJ takes the form as given in Eq. ( 2-28 ) with Zi mw given by 

fif2 r I'M 



(7-20) 
(7-21) 

(7-22) 



When Zi muJn are obtained, the energy and angular momentum fluxes averaged 
over t 3> At r are calculated by using Eqs. ( |2-31 ) and ( |2-32| ), respectively. Here, we 
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show these fluxes accurate to 0(e ) and to 0(v ) beyond Newtonian: 

(^E\ = E^+e 2 E^+0(e 3 ), (7-23) 
d -^) = j^ + e*jP + 0(e*). (7-24) 

We note that j] "* = E^/fl where Q = v 3 /M. We have already given the 5.5PN 
formula for E^ ' in section ||, Eq. (4-18). Hence our task here is to evaluate the 0(e 2 ) 
corrections. From the forms of r(t) and ip(t) given in Eqs. (7-5), we readily see that 
Zimuj n = 0(e'"') for n = ±1,±2. Thus, we only need the modes n = 0, ±1: As for 
the n = modes, we must retain terms up to 0(e 2 ), while for the n = ±1 modes, 
we only need terms up to O(e). Then the 4PN formulas for E^ and j {2) are found 
as 

e 2 E® =e 2 
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In Appendix |G|, we show each (£, m, n) component of the energy and angular mo- 
mentum fluxes. Note that there was an error in the coefficients of the e 2 v 4 terms in 



Ref. 17). This error is corrected in Eqs. (7-25) and (7-26) above. 

To express the energy and angular momentum fluxes in terms of the variable 
x = (M.!?^) 1 / 3 , we use Eq. (7-7). To 0(e 2 ), it can be easily solved for v as 
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(7-27) 



Then to 0(x 8 ) the energy and angular momentum fluxes are expressed as 
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where (d,J/dt) N is the Newtonian angular momentum flux expressed in terms of x, 

dJ A _32/M 2 Mi7| (7 .30) 



dt J N 5 \M 

and the e-independent terms in both (dE/dt) and {dJ /dt) are the same and are given 
by the terms in the case of circular orbit, Eq. (4-18), with the replacement v — > x. 

Finally, we consider the stability of circular orbits. We note that the following 
relation holds: 

§ ) - a, (f) - g - ^) , + art - + o( A 

(7-31) 

where u>± = mQ^^Qr and = Zt mw± . Here -H"(u) is an important quantity which 
determines the stability of circular orbits under the radiation reaction. Assuming 
the adiabatic evolution of the orbit, the evolution equations for tq and e due to the 
gravitational radiation reaction are written asIS 
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where 
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G(u) = £< 2 ) - Qjf) = H{v) - f{v)E^. (7-35) 
Using Eqs.(7-25) and (7-26), the 4PN formula of G(v) is calculated as 
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Note that [g(v)E^ + G{v)]/E^ 19/6 for w -> 0- i.e., in the Newtonian limit, the 
radiation reaction always reduces the eccentricity cfP . By a numerical calculation, 
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Apostolatos et al.cJ found that there exists a critical radius r c below which the 
circular orbit becomes unstable; r c ~ 6.6792M. On the other hand, we find the use 
of the 4PN formulas for E^ and G{v) gives r c ~ 7.38M. This indicates that a much 
higher order PN formula will be necessary to determine r c with good accuracy. 



§8. Slightly eccentric orbit around a rotating black hole 



In this section, we consider a slightly eccentric orbit on the equatorial plane 
of a Kerr black hole and calculate the leading order corrections of the eccentricity 
to the energy and angular momentum fluxes up to 0(v 5 ) beyond Newtonian. The 
calculation is parallel to the one given in the previous section. 

We consider the motion of a particle in the equatorial plane 6 = vr/2, hence we 
have C = 0. We define the radius r = ro as the one at which the potential R/r 4 is 
minimum; d(R/r 4 )/dr \ r=ro = 0. We define the eccentricity e such that r = ro(l + e) 
is a turning point of the radial motion at which R{r = ro(l + e)) =0. We assume 
e<l. Using these definitions of ro and e, E and l z are expressed as 

E = E® +eE^ + e 2 E^ +0(e 3 ), 
l z = lf) +d 0-) +e Hf)+O{e% 



where E^ and l z n ^ (n = 0, 1, 2) are given by 



1 - 2v z + qv 6 



(l-3v 2 + 2qv 3 )W 2 ) ' 
£«=0, 

{2) _ v 2 {\ - 3v 2 + qv 3 + q 2 v 4 )(-l + 6v 2 



+ 3q z v 



2(1 - 3v 2 + 2qv 3 ) 3 / 2 {-l + 2v 2 - q 2 v 4 ) 
rov(l — 2qv 3 + q 2 v 4 ) 



(1 - 3v 2 + 2qv 3 )( 1 / 2 ) ' 
^(2) Q r v 5 {q — 3v + qv 2 + q 2 v 3 ){— 1 + 6v 2 



2„,4\ 



3q z v 



2(1 - 3v 2 + 2qv 3 ) 3 / 2 (-l + 2v 2 - gV) 

where v = (M/ro) 1//2 . The post-Newtonian expansions of E^ and l{ n) up to the 
required order are 



M 3M 2 qM 5 / 2 
2r~ + ^f~^J 2 ~ 



' M 5M 2 3gM 5 / 2 \ 
2^"l^- + ^72-J 



+ o(A(s-i) 



l z = (Mro) 1 / 2 



1 + 



3M 3qM 3 / 2 t (27 
2r^ ZN 2 



+ Hr+<r -r 



M 2 15gM 5 / 2 



2r, 



5/2 



'g 2 M 2 3qM 5 / 2 \ 



2r 2 



2r, 



5/2 



+ o(^ 6 
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Now we solve the geodesic equations for a slightly eccentric orbit. The radial 
equation is 



We expand r(t) as 



Kit 



r(t) = r 1 + erW(t) + eV^(i) + 0{e 6 ) 



R_ 

J>2 ' 



•3) 



(8-4) 



and R/T 2 in terms of e and v using Eq. (8-1) and (8-2). Collecting terms of the 
equal order in e, we obtain 



^ 2 (1 " (r (1) ) 2 ), 



and 



1 drWdr( 2 ) 



r«f( 2 ) + go + gir w +? 2 (r w f, 



(i) 



^2 ^ 

where i? r , go, <7i and q<i are given in the post-Newtonian series forms as 



9o 



91 



92 



M 1 / 2 

3/2 



3M 3gM 3 / 2 (9 + 3g 2 )M 2 15gAf 5 / 2 6 

1 3/2 ol2 1 5/2 ^ I 



To 
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M 2gM 3 / 2 (6 + <? 2 )M 2 20gM 5 / 2 



3/2 



5/2 



+ 0(A 



2M 



2M 3oM 3 / 2 4M 2 60M 5 / 2 ^. f 
1 + + — -^ttt- + 0(^ e 



,3/2 



5/2 



3M 2<?M 3 / 2 (10 + g 2 )M 2 26gM 5 / 2 6 



(8-9) 
MO) 



We obtain r"(i) from Eq. ( |8-5| ) as 

r M(t) = cosi? r t, 



•11) 



where we set r(t = 0) = r (l + e). Substitution of Eq. <^TTj ) into Eq. (§]|) and 
yields after integration 



r (2) (t) = q 3 (l - cos Q r t) +94(1 -cos2J7 r t), 
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where q% = — qo and 54 = 52/2 . 

In the same way, we can solve the angular motion (p(t). From Eq. (2-18), we 
have dip/dt = <P/T, which can be expanded in terms of e using Eqs. (8-1), (8-2), 
), ( |8- 1 1| ) and ( |8-12j ). Integrating the resulting equation, we obtain 



cp(t) = Q(pt + epi sin Q r t + e 2 p2 sin Q r t + e 2 pz sin 2Q r t + 0(e 3 ), 



•13) 
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where 



4M 60M 3 / 2 (17 + g 2 )M 2 48gM 5 / 2 fiN 

Pi = -2- — + ^37^"- ^ + ^57^ + °^ )' 

'0 r ' '0 r ' 

2M 2qM 3 / 2 (1 - q 2 )M 2 6qM 5 / 2 fiN 
P2 = 2 + + ^ y + + 0(u 6 ), 



P3 



?'0 



„3/2 



„5/2 



6 M _ 2jM^ _ + HK-l^jg + 0( „ 6)i(8 . 14) 

4 4r r 3 ' 2 8r 2 8rg 



and 



3/2 



ro 



J 3 9M 9gM 3 / 2 3 (6 + g 2 ) M 2 60gM^ \ 6 







As in the case of the previous section, the fact that I? r 7^ JX, implies that these 
eccentric orbits are not closed. 

Using the above solution of the geodesic equations, we evaluate the source term 
of the Teukolsky equation. We set = tt/2 in the expressions of A nn Q etc. in 
Eqs. (2-25). Again, parallel to the discussion in section |7[ the orbits of our interest 
have the properties, 



r(t + At r ) = r(t), ^ 



t=t+At r 



dip 
~dt 



t=t 



where At r = 2ir/£2 r . Hence Eq. ( 2-26 ) reduces to the form, 

Zgnuo = Z£ mu] 5(lO —u> n ), 

where 

u> n = nf2 r + mf2 v , (n = 0, ±1, ±2, . . .), 

-At r 
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and 



Z, 



•19) 



ZiuJ -e lmuJ Jo 

with Wimu given by Eq. (2-27). 

Using the solution of the geodesic equation for r(t), we expand Wi mw in terms 
of e. The result takes the form, 



W, 



fo + ei /irW + h 



dt 



e 2 { / 3 r< 2 > + fj^ + / 5 (r«) 2 



dt 



' dr {i y 
dt 



+ h + °( e 
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where To ~ /s are time-independent coefficients. Inserting this form to Eq. ( |3-19 ) 
we obtain 



jjLTT 



iu n B mc 



£muj n 



2 2 

m pf 



imQ rPl Q 2 
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fomp 3 



+ /3<?3 



73 94 



+ ^hw + -f 2 mpiw - i f^q A w + ( 
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where <5„ )ri ' is the Kronecker delta. We see from this equation that Zi muJn = 0(e' n ') 
just as in the Schwarzschild case. Therefore we only need to retain the n = 0, ±1 
modes to evaluate the luminosity up to 0(e 2 ). 

We calculate the energy and angular momentum fluxes to 0(v 5 ) beyond the 
quadrupole formula and to 0(e 2 ) in the eccentricity. The time-averaged energy and 
angular momentum fluxes are given by Eqs. ( 2-31 ) and ( 2-32 ), respectively. In order 
to express the post-Newtonian corrections to the luminosity, we define r\i mn as 



d 4) -U d 4) , 



dt J 



Iran 



dt J 



N 



•22) 



where (dE/dt)^ is the Newtonian quadrupole luminosity given by Eq. ( 4T9 ). In the 
following, we show rj£ mn for m > modes, rjg^n for m < are obtained from the 
symmetry r\(,^ m ,n = Vt-m-n, which follows from the property of u> n in the present 
case, given by Eq. ( |7-21| ). 

For £ = 2, the 2.5PN formulas for rj£ mn are found to 



772,2,0 



107 v 2 
21 



+ 4ttv 3 - 6qv 3 + 



4784 n 2 4 

h 2 q 2 v 4 

1323 



428 irv 5 4216 qv 5 
+ 



21 



189 



*' As mentioned in the previous section, we have detected an error in the formula for 172,2,0 in 
Ref. 0). The term - 14270/ 147e 2 v there is in correct. Accordingly, formulas for dE/dt and dj z /dt 
below are also corrected in this paper. 
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+e 2 - 10 + ^ 46ttu 3 + 84gv 3 — 23<?V 

4748 ttv 5 2675 qv 5 \ 
+ 21 + 189 ) ' 

2 /729 3645 v 2 2187 vrv 3 3645 gu 3 
772,2,1 ~ 6 ^"6T 64" + 32 32 

24057 v 4 2187 q 2 v 4 6561 ^v 5 9477gt> 5 \ 
+ — + ^ — + 



r?2,2 - 1 " ^ ^64 + ^48~ + 32 32 

2224681 v 4 99 q 2 v 4 615 vrt; 5 _ 27857 gt; 5 ' 
+ 112896 + 64 + 112 336 ^ 

_ v 2 qv 3 I7v 4 q 2 v 4 ttv 5 793qv 5 
V2,l,0 ~ oa io cn/i 1« 



256 64 16 112 

9 1041 v 2 9 ttv 3 153 qv 3 



36 12 504 16 18 9072 

2 ( -2v 2 2qv 3 93 v 4 q 2 v 4 19 ttv 5 27113g?; 5 
+e ^ 9 + 3 + 7l2 2 36 18144 

2 /4v 2 4gt; 3 172t; 4 2 4 2794gv 5 \ 
m,i,i = e — ^ u +^7— + 



r?2,0,±1 6 ^96 672 + 48 + 16 + 169344 32 

83vrt> 5 1255 qv 5 \ 
168 504 ) ' 

and 772,1,-1 becomes 0(v 6 ). Putting together the above results, we obtain (dE/dt)^ 
J2 mn {dE/dt) emn for I = 2 as 

dE\ fdE\ f \277v 2 3 73qv 3 37915v 4 

~dt) 2 ~ \~dt) N \ 252 + nV 12~ + 10584 

33 g 2 v 4 _ 2561 ttv 5 201575 qv 5 
+ 16 126 + 9072 

2 /37 2581 1; 2 1087 ttv 3 211 qv 3 325393 1; 4 105 q 2 v 4 
+e \24 252 + 48 6 + 21168 + 8 

29857 7r v 5 11293 gw 5 \l 
5i— *— «H ' (8 ' 23) 



9 3 63 " 9 567 
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For I = 3, we obtain 

1215 v 2 1215 v 4 3645 tt v 5 1215 qv 5 
773,3,0 ~ 896 112 + 448 112 



38 



2 / -10935 v 2 37665 v 4 142155 irv 5 134865 q v 5 
+e y 448 + 256 896 + 448 
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o v qv \ 
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2688 1008 



and 773,1,-1 becomes 0(v 6 ). Thus we obtain 



gLE\ _ fdE\ [ l367t; 2 32567 v 4 / 16403 tt 896 g \ 5 
~dt) 3 ~ \ dt) N \ 1008 3024 + V 2016 81/ " 
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+e ^ 252 864 + ^ 448 56 ) ) } A * > 



For t = 4, we have 
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e 2 v 4 



??4 ' ' ±1 " 225792 ' 
and 774,2,1 becomes 0(v e ). Hence we have 



dE\ _ fdE\ J 8965 t; 4 2946739 e 2 v 4 



dt J 4 \dtJ N \ 3969 127008 J 1 y 

Finally, gathering all the above results, we have the luminosity up to 0(v 5 ) as 

dE\_(dE\ L 1247 v 2 3 73gt> 3 44711 v 4 33 g 2 1; 4 

di/ W*/n\ 336 + nV 12 9072 H 16 - 

8191 vru 5 3749 qv 5 2 / 37 65 u 2 1087 irv 3 211 qv 3 

672 h 336 +g ^24 21~ + 48 6 

465337 v 4 105 g 2 v 4 118607 ttv 5 95663 gv 5 \l 

9072 + 8 1344 672 J J ' ^ 8 ' 26 ^ 

If we set q = 0, the e 2 correction terms in the above formula completely agree with 
the corresponding terms in Eq. (7-25) in the previous section. 

To compare our results with those derived in the standard post-Newtonian 
method, it is convenient to change the parameter from v to x = (MQ^) 1 / 3 . The 
relation between v and x is given by 



1 + -x 3 + e 2 (- - -x 2 + -ox 3 - 6x 4 - q 2 x 4 + — qx 5 
3 1 2 2 3 y H 2 H 



1-27) 



Then we obtain 



dE\ dE\ f 1247 „ 11 „ „ 44711 4 33 . „ 

( — ) = — < 1 x z qx* + 4ar i 7r x H x 4 <T 

\dt/ \dt J N \ 336 4 9072 16 

59 r 8191 ,5 /157 6781 2 2009 , 2335 , 

qX X 7T + X QX H X 7T 

16 672 V 24 168 72 4 48 

14929 . 281 4 „ 2399 = 773 «, \ 9 1 , 

-Tsr^+ii-^-s-^-r^}- (8 ' 28) 

where (dE/dt) N is the quadrupole flux expressed in terms of x, Eq. (|6-6|). We find 
that the termswhich are proportional to e 2 agree with the formulas derived by Peters 
and Mathews Ea) at leading order, Galt'sov et al.Q) and Blanchet and Schafer at v 2 
order|, Blanchet and Schafer at v 3 order for q = 00 and Shibata at v 3 order for 
g 7^ 0E3*, if we expand their formulas by e assuming eCl and /i/M <C 1. 

From Eq. ( 2-32j ) , the partial mode contributions to the angular momentum fluxes 



for I = 2, 3 and 4 are calculated to be 

fdJ z \ /dJ z \ ( 1277 v 2 3 61 q v 3 37915v 4 33q 2 v 4 
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where (dJ z /dt)^ is defined by 
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Total angular momentum luminosity is then given by 
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The e 2 terms in the above also agree with the corresponding terms in Eq. (7-26). 
The angular momentum flux expressed in terms of x = (MQ^) 1 / 3 is given by 



dJ z 
dt 
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dt 
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where (dE/dt)^ is the Newtonian flux expressed in terms of x, Eq. (7-30). 

§9. Circular orbit with small inclination from the equatorial plane 



In this section, we consider the case of a circular orbit at r = ro with small 
inclination from the equatorial plane. We evaluate (dE/dt) and (d.J z /dt) to 0(v 5 ) 
beyond Newtonian. 
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In this case, the orbital plane precesses around the symmetric axis. The degree 
of precession is determined by the value of the Carter constant C. If ro and C are 
given, the energy E and the z-component of the angular momentum l z are obtained 
by the two equations, R = and dR/dr = 0, where R is a function defined by 
Eq. ( 2TSp . We introduce a dimensionless parameter y defined by 



y 



-■ Q 2 

Q 2 ' y 



H + a 2 (l 



E 



(9-1) 



We assume y is a small number. Since Q 2 ~ I 2 and C ~ I 2 + Z 2 , this is physically 
equivalent to assuming l 2 + l 2 <^ I 2 . Since we do not need the exact expressions for 
E and l z in terms of ro and y, we show them to the first order in y as well as to 
0(v 5 ). They are given by 



E = 1 



M 3M 2 
2r~^ + lkf 
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l5M 3 / 2 a 
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2r, 



(l-y)+0(v e 



(9-3) 



where note that a = Mq (\q\ < 1). 

To solve the geodesic equations under the assumption y <C 1, we first set 6 
tt/2 + y l l 2 9' and consider the geodesic equation for 9. It then becomes 



dr J 



E 2 



Q 1 



Sill' 



y 



a 2 (l 



E 2 ) + 



I 2 



cos 2 {y 1 / 2 e > ) 



(9-4) 



Since the right hand side of Eq. ( |9-4j ) contains only even-functions of y 1//2 0', we can 
solve it iteratively by expanding 9' as 

0' = 9 {0) + y9 {1) + y 2 9 {2) + ■■■. (9-5) 

This method is similar to the one we have used in section [7| or ^. However, here we 
only consider the lowest order solution 9(q\- This means we take into account the 
effect of inclination up to 0(y), as seen from the structure of the geodesic equations. 
The equation for 0( O ) is 



( d9 {0) 
V dr 



or dividing it by (dt/dr) 



d9 



(0) 



dt 



Q 2 



0~ 



7 (o) 



where 



a = —a(aE — l z ) + 



A(r ) 



rl + a 2 ) 



alz}- 



(9-6) 



(9-7) 



(9-8) 
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Then the solution is easily obtained as 



ho) = sin(^); Q e = — , 



where we have chosen 9 



(o) 



at t = 0. Thus we have 
0= J +y 1/2 sin(^). 



(9-9) 



(9-10) 



Note that the solution ( 9-10D implies that the inclination angle 9{ is indeed given by 
9i = y 1 / 2 in the present approximation. 

Next, we consider the geodesic equation for cp. Taking account of the terms up 
to 0(y), it becomes 



dtp k 
dt a 
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The solution to Eq. (9T1) with <p = at t = is 

ip = QJ, - 2/777- sin(2 fl e t). 



(9-11) 
(9-12) 

(9-13) 
(9-14) 



Note that JX> ^ ^e- This means the precession of a test particle orbit around the 
spin axis of the black hole. Specifically, to the order required for the present purpose, 
we have 





M 1 / 2 




3/2 






(2 e = 




3/2 
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M 1 / 2 a 3 (M l l 2 a a 2 



r 3/2 + 2 V 




V r 3 J 2 
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+ 0(v 6 ) 



' 3M x / 2 a 3a 2 6 , ^. . 



(9-15) 



We see that £2^ — Qq — ► 2Ma/ri for tq —>■ oo and y — > 0, which is just the Lense- 



Thirring precessional frequency!-^. 

Now we are ready to calculate the source integral for the amplitude Ze muJ . Anal- 
ogous to the case of an eccentric orbit considered in section [7| or ||, Eq. (2-26) can 
be simplified further by noting that the orbits of our interest have the properties, 



9(t + At e ) = 9(t), 



<p(t + At e ) = <p(t) + A<p, 



(9-16) 
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where Atg is the orbital period of the motion in the ^-direction and Aip is the phase 
advancement during Atg. In other words, we have 

27T ,.. Aip 



Then we obtain 

where 
and 



Z 



£mu)„ 



Atg' " V Atg " 

(n = 0,±1,±2 
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nQg + mfi v 



£moj n j 



(9-17) 

(9-18) 

(9-19) 
(9-20) 



2^n^„ JO 

with W tmuJn being given by Eq. (2-27). 

Let us discuss the final form of Zi muJn . In the present case, up to 0(y) the 
integrand We muin has the form, 



W, 



i = go + y 1,2 giO{o) + V920f 0) 

d6/n\ _ d9( n \ f d6(Q) x 2 



y /2 ^ + ^ {0 )^ + ^^) + 0{y*l% (9-21) 



where go ~ g§ are complicated functions of rg. Using an approximation, 
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we have 



At 
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Thus the amplitude ^£ mWn is found to have the form, 



J £mui n 



(9-23) 
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(9-24) 



where Z l > 3 are functions of rg. Here, it is worth noting the symmetry of Z^ muJn . The 



spin weighted spheroidal harmonics have a property _2S^ n (6>) 



(-iy-2sr m n (n- 



Then, from Eqs. (Q) and (^), we have Z £ _ mw _ n = (-l) n+£ Z ta 
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Now we evaluate the energy and angular momentum fluxes at infinity. The en- 
ergy and angular momentum fluxes averaged over t ^> Ate are give by Eqs. ( |2-31| ) 
and ( |2-32| ), respectively. Then we see from Eq.(9-24) that the n = ±2 modes con- 
tribute to the luminosity at 0(y 2 ). Thus, when we calculate the luminosity to 0(y), 
we need to include only the n = 0, ±1 modes. In order to express the post-Newtonian 
corrections to the luminosity, we define rj£ mn as 

dE\ 1 fdE\ 



dE \ ~ 1 ( dE \ T) (9 25) 

dt J Iran 2 \ dt J w 



where (dE/dt)^ is the Newtonian quadrupole luminosity given by Eq. Q4-19 ). 

For £ = 2, the results are as follows. If \m + n\ >2orm + n = 0, % mn becomes 
of 0(v e ) or higher. The remaining r\i mn which contribute to the 2.5PN luminosity 
formula are given by 
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21 1 1323 2 
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v b H 
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Putting together the above results, we obtain (dE/dt)i = J2mn(dE/dt)g mn ioi £ = 2 
as 



12 Oil = V ( - l:qv 3 - ^ + + ^ - ^9» 5 I ■ P'M) 



diy 2 V^/n\ 252 12 y V 2) 



37915 4 
10584 W 
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H g V gVy v 5 H gu 5 1 - - H . (9-27) 

16 96 126 9072 y V 2/ ' 

For £ = 3, the non-trivial r/£ mn are given by 

1215 2 1215 4 3645vr = 1215 = 

fh+3 n = v v H v qv 

1 896 112 448 112 H 

( 3645 2 3645 4 10935tt = 3645 t ' 
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tf V 1792 224 896 112 y 

5 4 

%±3=fi = -^v y, 



Chapter 1 Black Hole Perturbation 



45 



V3±2 
V3±2±l = y 

%±2=fi = y 
m±i o 



5 4 40 . / 20 4 100 

— v qv + y\ v -\ qv 

63 189 y tf V 63 189 y 
Z3645 2 3645 4 10935vr = 6075 



1792 224 

5 2 5 

v 

16128 



896 

4 , 57r 
-f + 



qv 



v 5 + 



224 

25 «= 



1 



8064 

+y 



V 



3024 8064 18144 

1 4 7T c 17 e 



11 



1512 



4032 



16128 



11 

3024 1 



9072 



llvr 
8064 



v° + 



95 
9072 



qv 



m±i±i = y\ 
m o±i = y 



f25 4 

v 

V126 

1 , 



qv 



2688 



80 
189 

504 1344 



7T 



11 

1008 



The other r]£ mn are of 0(v G ) or higher. Then we obtain 
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For £ = 4, we have 



??4±4 o 
??4±3±1 
??4±3t1 
??4±2 



1280 

567 
2560 

"567 



v\l - 2y), 



v y, 



1134 
5 



v y, 



3969 



v\l-8y), 



v y, 



and the others are of 0(v G ) or higher. Hence we obtain 
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Finally, gathering all the terms, the total energy flux up to 0(v 5 ) is found to be 
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96 * " 672 336 

Using the above results for % mn , the time-averaged angular momentum flux is 
calculated from Eq. ( 2-32 ). The partial mode contributions of the I = 2, 3 and 4 
modes are calculated to give 
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where (dJ z /dt)^s is defined in Eq. ( |8-29D . The total angular momentum flux is then 
given by 



dJ z 
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(9-34) 



We note that the result is proportional to (1 — y/2) in the limit q — > 0. This is 
simply because the orbital plane is slightly tilted from the equatorial plane by an 
angle 9i ~ y 1 ^ 2 , hence dJ z /dt ~ (d Jtot /dt) cos 0j . 

§10. Adiabatic backreaction 



In the preceding sections, we have evaluated the energy flux (dE/dt) and the 
z-component of the angular momentum flux (dJ z /dt) emitted to infinity by a particle 
for various cases. By emitting gravitational waves, a particle orbit will suffer from 
radiation reaction. In the limit of small fJ>/M, the reaction time scale will be much 
longer than the characteristic orbital time scale; t reac t ~ M 2 /fi S> At. Hence the 
evolution of the orbit will be well described by the adiabatic backreaction. 

In the case of orbits around a Schwarzschild black hole or orbits confined on 
the equatorial plane around a Kerr black hole, it is straightforward to calculate 
the evolutionary path under radiation reaction because the orbits are completely 
specified by the energy E and the z-component of the angular momentum l z , hence 
their time derivatives can be simply evaluated by equating them with —(dE/dt) 
and —(dJ z /dt), respectively. However, once we consider motions off the equatorial 
plane of a Kerr black hole, the orbits cannot be specified by E and l z alone but the 
specification of the Carter constant C becomes necessary. Unlike E or l z , since C 
is not associated with the Killing vector of the spacetime, one cannot calculate the 
radiation reaction to C by simply calculating the gravitational waves at infinity. This 
implies that we have to derive a local radiation reaction force term to the geodesic 
equation by evaluating the metric perturbations around the particle, as is done in the 
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derivation of radiation reaction force in the standard post-Newtonian method. For 
almost Newtonian orbits, applying a post-Newtonian radiation reaction force, Ryan 
derived the evolution equation for the Carter constant c3> . However, no relativistic 
treatment has been done so far. This is a challenging issue. An approach to this 
issue is discussed in Chapter 7. 

In this section, instead of attacking this very difficult problem, we discuss some 
general properties of the adiabatic radiation reaction in a restricted class of orbits. 
Namely we consider orbits which are circular or those having small eccentricity. We 
clarify the conditions for circular orbits to remain circular under radiation reaction. 
A detailed discussion on this matter has been given by Kennefick and OriE^. We 
give a less detailed but more general discussion below. 

We recall that the radial velocity u r = dr/dr is written in terms of the first 
integrals of motion in the test particle limit as 

{Eu r f = R{I\r\ (10-1) 

where P = (E,l z ,C) and R(P,r) is independent of 8 and <j>. First let us consider 
orbits which are circular in the test particle limit. These orbits are determined by 
the conditions, 

f) P 

R(P,r) = 0, ^-(P,r) = 0. (10-2) 

Eliminating r from these equations gives an implicit relation among P's. For exam- 
ple, 

f(P) = R(P,r(P)) = 0, (10-3) 
where r(P) is obtained by solving the second of Eq. (10-2) for r. This equation 



determines a two-dimensional hypersurface S in the 3-dimensional space M of P's. 
The adiabatic evolution of an orbit is characterized by slow evolution of P, i.e., 
P = 0((a), where /i is the mass of the particle. Then a necessary condition for 
circular orbits to remain circular under radiation reaction is that we have f(P) = 
fiP = 0(fi 2 ). In other words, the vector P on S is tangent to S to O(fi). This 
condition can be shown to hold by the following theorem. 

Theorem: If the radiation reaction to the r-component of the acceleration is of order 
fi; a r := du r ' jdr = O(fi), i.e., the r-component of the radiation reaction force is 
well-defined and finite, then for orbits which are circular in the test particle limit; 
i.e., u r = O(n), the radiation reaction to P is constrained by the equation, 

f) P 

Wi (P,r)P = 0(n 2 ), (10-4) 

where the argument r is to be replaced by r(P) after differentiation. 

Proof: It is almost trivial. Just taking the T-derivatives of both hand sides of 
Eq. (|iol ) gives Eq. (|l0l ). Q.E.D. 



Thus, since 

• / dR OR ■ dr\ ■ 

fin = ( W P,rF)) + ^ P ^n) W ) I\ (10-5) 
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and the second term in the parentheses vanishes by definition, we have / = 0(fj, 2 ). 

This theorem alone, however, does not mean that circular orbits remain circular, 
since we have constrained the first integrals to be those for circular orbits from the 
beginning. Let us explain the reason. Since we may regard P a vector field in M, 
what we need for circular orbits to remain circular is the regularity of P in the vicinity 
of the hypersurface S. In other words, if the vector field P is not differentiable on S, 
an orbit on S may spontaneously deviates away from S. A simple illustrative example 
is the case I± = vXl at ij_ = where I± is the component of P perpendicular to S. 

Thus, provided P is regular in an open neighborhood of S, the above theorem 
implies that a circular orbit in the test particle limit remains circular under adiabatic 
radiation reaction. In this case, the radiation reaction to the Carter constant, C, is 
determined by the radiation reaction to the energy, E, and the z-component of the 
angular momentum, l z . Specifically we have 

C = (-^ (E (r 2 + a 2 ) - al z ) (r 2 + a 2 ) - 2a (Ea - l t )\ E 

+ (-2-| (E (r 2 + a 2 ) - al^ +2{Ea- l g )j l z . (10-6) 

Yet this is not the end of the story. What we have shown is that P lies on S. 
But if P slightly off the hypersurface S is diverging away from S, circular orbits 
will be unstable. Thus the condition for the stability of circular orbits is that S 
is an attractor plane of the vector field P. However, the notion of divergence or 
convergence of a vector depends on the metric of the space M, but we have no guiding 
principle to determine the metric. This implies that the notion of the attractor or 
the stability is ambiguous. 

Nevertheless, extrapolating from the case of Newtonian orbits, there seems to 
exist a natural choice of the metric. Namely, as the distance of the orbit from the 
hypersurface S of circular orbits, we define the eccentricity of an orbit as given 
in sections |?] and ||. With this choice of the metric, let us consider the adiabatic 
radiation reaction problem in more specific terms. 

Let us parametrize an orbit in terms of the mean radius ro, the eccentricity e 
and the square root of the Carter constant y := C 1 / 2 , instead of the energy E, the 
angular momentum l z and the Carter constant C . The mean radius ro is defined by 
the equation, 

= R'(I\r ), (10-7) 

where the prime denotes the partial derivative with respect to r. This definition says 
that r is maximum at r = ro. The eccentricity e is defined by setting the maximum 
radius to r = ro(l + e), i.e., 

= R(P,r (l + e)), (10-8) 

This definition guarantees that e = corresponds to a circular orbit. Assuming 
e<l, the above equation can be expanded in powers of e as 

= R(P,r ) + l -R"{P,r Q ){r e) 2 + ±R®{r,r )(r e) 3 + j { R {i \l\ r )(r e) 4 + • • • , 

(10-9) 
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where R^ is the n-th derivative of R with respect to r. The parameters (ro,e, y) 
are chosen because the geodesic trajectory x M = z^(t) allows perturbative expansion 
in powers of e and y at least for eCl and y <C 1. Therefore the first integrals of 
motion P = (E,l z ,C) will be regular functions of the parameters (ro,e,y). On the 
other hand, if we consider (ro, e, y) as functions of /*, it should be noted that e is not 
a regular function of P in a neighborhood of circular orbits because of the absence 
of a term linear in e in the right hand side of Eq. ( |10-9| ) . 

Now we consider the adiabatic evolution of e under the radiation reaction. Tak- 
ing the r-derivative of Eqs. ( 10-7 ) and ( 10-9j ), we get 

= R' ,P + fl&Vo , 







R^P + flg (r r e 2 + r 2 Q ee) + X - (r^P + R^ro) r 2 e 2 
+\R i i ) rle 2 e + ±i#>rg C 3 e + O (e 3 , ee 4 ) , 
where -R 0i = d 2 R/dPdr etc. Equation (10T0) determines r'o as 



o,; 



P. 



Substituting this into Eq. (10-11), we obtain the expression for e as 



" pr 2 nil 

er n 



+e z 



Ro,i + 



roX,i+ko+ 2 



e rp_Rp_ 
2 itf 



(3) 



0,t 



2 o(3) ~ 



li?, 



(4) 



1 / R, 



(3)' 



6 i?" 



i? ' 



(10-10) 



(10-11) 



(10-12) 
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(10-13) 



Since the trajectory z^{t) is assumed to be analytic in (ro, e, y), it is reasonable 
to further assume that P are regular functions of (ro, e, y). Then we can expand P 
with respect to e as 

P (r , e, y) = P^(r , y) + eZ^ro, y) + e 2 P^(r , y) + • • • . (10-14) 

Then we obtain 
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+0(e 2 



(10-15) 



As shown by the theorem above, Eq. ( 10-4 ), the leading term of order e 1 vanishes 
provided the radiation reaction force is finite: 



R,J m = 0. 



(10-16) 



The next order term determines whether the circular orbit remains circular or not. 
If it does not vanish, the eccentricity will spontaneously develop as 



(10-17) 



Here the regularity of P comes into play. As noted above, e is singular on the 
hypersurface S. Hence if P is regular on S, /W(ro,y) should vanish. By a detailed 
analysis, it is shown in Ref. ^) that this is indeed the case. The physical reason is 
rather simple: If one considers a slightly eccentric orbit, there appears a frequency 
of wobbling motion due to the eccentricity, say f2 e . In general the ratio of Q e to the 
frequency of the motion in the 9 or ip direction is an irrational number. Hence the 
part of the metric perturbation which is proportional to e will have frequencies that 
are integer multiples of J? e , and the same property is shared by the corresponding 
term of the backreaction force linear in e. Since any sinusoidal oscillation has zero 
mean when averaged over time longer than its period, this implies there will be no 
term linear in e in the adiabatic expression of P. 
Thus we have 



+ [r + 



\_rlBf 
2 R o , 



r> m 



e + 0(e 2 ), (10-18) 



and circular orbits will remain circular under radiation reaction. As for the stability 
of circular orbits, whether the eccentricity decreases or increases is determined by the 
sign of the coefficient of e in the right hand side. Thus it is necessary to calculate the 
radiation reaction to the Carter constant to determine the stability. As mentioned 
in the beginning of this section, this is a challenging issue. Finally, we should again 
note that the meaning of stability does depend on the definition of the eccentricity, 
i.e., how we define the distance from the hypersurface of circular orbits S. 



§11. Spinning particle 



So far we have considered only a monopole particle orbiting a black hole. How- 
ever, in a realistic binary system of compact bodies such as a neutron star-neutron 
star, black hole-neutron star or black hole-black hole binary, both bodies may have 
non-negligible spin angular momenta. Hence it is desirable to take into account not 
only the spin of a black hole but also the spin of a particle in the calculations of 
gravitational waves from a particle orbiting a black hole. 

To incorporate the spin of a particle, one must know (1) the equations of motion 
and (2) the energy momentum tensor of a spinning particle. Fortunately, we know 
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that (1) have been derived by Papapetroulll' , Dixon© and Wald© and (2) has also 
been derived by Dixon c3> . Hence, by using the expression for the energy momentum 
tensor of a spinning particle as the source term in the Teukolsky formalism 0>, we can 
calculate the gravitational waves emitted by a spinning small mass particle orbiting 
a rotating black hole. One may regard this particle as a model of a small Kerr 
black hole, but it may be appropriate here to give a word of caution. A Kerr black 
hole of mass (i and the spin parameter S, where S is defined so that fj,S gives the 
spin angular momentum, has quadrupole (I = 2) and higher multipole moments 
(£ > 2) proportional to fiS e as well. Since we neglect the contributions of these 
higher multipole moments here, our treatment will be valid only up to O(S) if we 
regard the particle as a Kerr black hole. To incorporate the contributions of all 
higher multipole moments to represent the Kerr black hole is a future problem to be 
investigated. 

Here we review the results obtained by Tanaka et al. i3. We concentrate on 
the leading effect due to the spin of the small mass particle. We consider a class of 
circular orbits which stay near the equatorial plane with the inclination solely due to 
the spin of the particle, i.e., those orbits which would be confined in the equatorial 
plane if the spin were zero. Then we calculate the gravitational wave luminosity to 
0(f 5 ) with linear corrections due to the spin. 

11.1. Equation of Motion and Source Term of a spinning particle 

To give the source term of the Teukolsky equation, we need to solve the equa- 
tions of motion of a spinning particle and also to give an expression for the en- 
ergy momentum tensor. In this section we give the necessary expressions, following 



Refs. |g),|2D 



Neglecting the effect of the higher multipole moments, the equations of motion 
of a spinning particle are given by 



^t) = --R» vp Mt))v v {t)S<>°{t) 



^-S^{t)=2 P ^{t)v u \t), (11-1) 
dr 

where v^{t) = dz^(r) /dr, r is a parameter which is not necessarily the proper time of 
the particle, and, as we will see later, the vector p^(r) and the antisymmetric tensor 
S^ u (t) represent the linear and spin angular momenta of the particle, respectively. 
Here D/dr denotes the covariant derivative along the particle trajectory. 

We do not have the evolution equation for u**(t) yet. In order to determine 
v m (t), we need to i mp ose a supplementary condition which determines the center of 
mass of the particlelH2P , 

S^(t) Pu (t) = 0. (11-2) 

Then one ran show that p^ = const, and S^ U S^ = const, along the particle 
trajectory©. Therefore we may set 

= put* , u ^ = -1, 
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S 2 = S ^ = —S^S*", (11-3) 



1 

27? 

where fi is the mass of the particle, is the specific linear momentum, and is 
the specific spin vector with S its magnitude. Note that if we use instead of S^" 
in the equations of motion, the center of mass condition ( 11-2| ) will be replaced by 
the condition 

PnS» = 0. (11-4) 

Since the above equations of motion are invariant under reparametrization of 
the orbital parameter r, we can fix r to satisfy 

""(tKCt) = -1. (11-5) 

Then, from Eqs. (11-1), ( ll-2j ) and ( |11-5| ), v^{t) is determined asS* 

v »(t)-u»{t) = ~ (m 2 + \r x(Cv (z(t))S x Ht)S^(t^ 1 S^{t)R V(XTK (t)uO(t)S™{t). 

' (H-6) 
With this equation, the equations of motion (11-1) completely determine the evo- 
lution of the orbit and the spin. Note that = + 0(S 2 ), hence and are 
identical to each other to O(S). 

As for the energy momentum tensor, Dixon gives it in terms of the Dirac 
delta- function on the tangent space at x^ = z^(t). For later convenience, in this 
paper we use an equivalent but alternative form of the energy momentum tensor, 
given in terms of the Dirac delta-function on the coordinate space: 



T a ?(x) =J drl p( a (x, t)v® (x, t) Si4) 



-V 7 ^(,,rK)(x,r) j(4)( ^ (T)) )|, (11-7) 



where v a {x,r), p a {x : r) and S a ^(x,r) are bi-tensors which are spacetime extensions 
of v^(t), P^(t) and S^ u (t) which are defined only along the world line, x^ = z^(r)0 
To define v a (x, z(t)), p a (x, z(t)) and S a ^(x, z(t)) we introduce a bi-tensor g a ^(x, z) 
which satisfies 

hmg^(x,z(r))=5 a IM , 

KmV f3 g a ^x,z(T)) = 0. (11-8) 

For the present purpose, further specification of g a ^{x, z) is not necessary. Using this 
bi-tensor g a ^(x,z), we define p a (x, r), v a (x,r) and S a ^(x,r) as 



p [x,r = g 



v (x,z{T)y(T), 



*' In the rest of this section, we use fi, u, <7, • • ■ as the tensor indices associated with the world 
line z(r) and a, ft, 7, • • ■ as those with a field point x, and suppress the coordinate indices of z(r) 
and x for notational simplicity. 
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S°0(x, r) = g% (x, z{t)) f u (x, z(r)) S^(r). 



(11-9) 



It is easy to see that the divergence free condition of this energy momentum 
tensor gives the equations of motion (11-1). Noting the relations, 



V f5 g a ,(.x,z(T))5^(x,z(r)) = 0, 



v a (x)V c 



dr 



(11-10) 



the divergence of Eq. (11-7) becomes 

r / s 8^(x — z(t)) / fj i , . 

V fi T*(x) = J drg a ,(x,z(r)) ^ J + ^R^(r))v u (r)S^(r) 

8^ (x - z{t)) 



+- dTvAg^{x,z{r))f v (x,z{r))- 



x (—S^(t)-2 P ^{t)v v \t) 
\cIt 



fll-111 



Since the first and second terms in the right-hand side must vanish separately, we 
obtain the equations of motion (11-1). 

In order to clarify the meaning of p^ and S^ v , we consider the volume integral of 

this energy momentum tensor such as / gJ l T al3 dU/3, where we take the surface 

S{tq) to be perpendicular to u a (ro). It is convenient to introduce a scalar function 
t(x) which determines the surface 2J(tq) by the equation t(x) = tq, and dr/dx 13 = 
—up at x = z(tq). Then we have 

^ gfT^dEf, = J d A x^-g^5{r{x) - r )g a ^{x) 



J dr'l 8{t' - t ) 



1 Du u 
2~dr~ 



^(ro). 



(11-12) 



where we used the center of mass condition and the equation of motion for S^ v . We 
clearly see p^ indeed represents the linear momentum of the particle. 

In order to clarify the meaning of S^ v , following Dixon W, we introduce the 
relative position vector 

X" :=-g^d v a(x,z), (11-13) 

where a(x, z) is the squared geodetic interval between z and x defined by using the 
parametric form of a geodesic y(u) joining z = y(0) and x = y(l) as 

dy a dyP 



a[x, z) 



2 Jo ^ al3 du du 



■du. 



(11-14) 
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Then noting the relations 

limX^ = 0, limX fi R = 5 l i, (11-15) 

X— >z x—>z >" P 

it is easy to see that 

= 2 f X^g^T^dE*. (11-16) 
Je to 

Now that the meaning of S^ u is manifest. Prom the above equation, it is also easy 
to see that the center of mass condition ( jll-2 ) is the generalization of the Newtonian 
counter part, 

' d 3 xp(x)x i = 0, (11-17) 



where p is the matter density. 

Before closing this subsection, we mention several conserved quantities of the 
present system. We have already noted that p^ = —p? and S^S^ = S 2 are 
constant along the particle trajectory on an arbitrary spacetime. There will be an 
additional conserved quantity if the spacetime admits a Killing vector field 

=0- (H-18) 

Namely, the quantity 

Qr-=p^-^^, (ii-i9) 

is conserved along the particle trajectory^). It is easy to verify that is conserved 
by directly using the equations of motion. 

11.2. Circular orbits near the equatorial plane 

Let us consider circular orbits around a Kerr black hole with a fixed Boyer- 
Lindquist radial coordinate, r = r . We consider a class of orbits that would stay 
on the equatorial plane if the particle were spinless. Hence we assume that 8 := 
9 — 7r/2 ~ 0(S/M) -C 1. Under this assumption, we write down the equations of 
motion and solve them up to the linear order in S. 

In order to find a solution representing a circular orbit, it is convenient to intro- 
duce the tetrad frame defined by 

*\ = 

e%= ( — =hhj 6*. 0. U. - — —>.\u(f). (.11 -20) 

where E = r 2 + a 2 cos 2 6, and e a ^ = (e a t , e a r , e a e , e a v ) for a = ~ 3. Hereafter, we 
use the Latin letters to denote the tetrad indices 




Chapter 1 Black Hole Perturbation 



55 



For convenience, we introduce uj\ ~ ujq to represent the tetrad components of 
the spin coefficients, u) ab c = e a fJ, e b l/ e c U] ^, near the equatorial plane: 

n l ^/n2\ a2 — Mr 

w oi = w oo =ui + 0(d ), OJi 



w 31 — w 30 — w 13 — w 10 — w 03 



r 2^1/2 
3 



u 01 6 = uj 2 + O(0 z ), uj 2 := -j 



1 _ 2 _ 1 

^22 — _w 21 — ^33 



_ 2 _ 1 
w 02 — w 00 — w 12 



2 ) 
Z\V2 



-u 31 6 = u 3 + O(r), w 3 := — 



-a; 



ii 



w 4 + 0(6> 2 ), cj 4 : = 



a 



w 32 — w 30 







w 23 

w 33 2 = -w 32 3 = flw 6 + O(0 2 



u 20 — w 03 — _w 02 



to 5 + O(0 2 ), ui 5 :~- 



CJ 6 



(r 2 + a 2 ) 



(11-21) 



Since the following relation holds for an arbitrary vector 

»dr J dr 1 "bcVJ, 
the tetrad components of Df^/dr along a circular orbit are given explicitly as 



e° — P 


= /°- 


(^/i + ^C/ 2 ) +O(0 2 ), 




»dr J 


= / 1 - 


(^/° + J B/ 3 + i?/ 2 ) +O(0 2 ), 




e 2 


= / 2 - 


(flC/ + #£>/ 3 - S/ 1 ) + 0(# 2 ), 




e 3 — P 


= / 3 - 


(-Bf-eDf^+oie 2 ), 


(11-22) 



where A, B, C, D and E are defined by0 



A 


:= lo iV ° 


+ 002V 3 


B 


:= UJ2V 


+ UJ3V 3 


C 


:= 


+ u; 5 v 3 


D 


■= U5V 




E 


■= w 3 v 2 





(11-23) 

and we have assumed that v 1 = and v 2 = O(0). 

For convenience, we rewrite the equations of motion by changing the spin vari- 
able. Instead of the spin tensor, we introduce a unit vector parallel to the spin, ( a , 
defined by 

C := ? = -^e a fec yS crf , (11-24) 



S 2/iS 



*' The symbols A ~ _B used here to define the auxiliary variables are applicable only in this 
subsection, and not to be confused with quantities defined with the same symbols such as E for 
energy, in the other sections. 
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or equivalently by 

S ab = fiSe ab cd u c ( d , (11-25) 

where e abcd is the completely antisymmetric symbol with the sign convention £0123 = 
1. As noted in the previous subsection, if we use the spin vector as an independent 
variable, the center of mass condition (11-2) is replaced by Eq. ( 11- 4j ) , that is 

( a u a = 0. (11-26) 

Then the equations of motion reduce to 

rlii a 

= co bc a v b u c - SR a , 

dr 

d( a 



where 



<h Lu bc a v b ( c - Su a ( b R b , (11-27) 



R a := R*\ cd v b u% d = ^R" bcdV b S cd , (11-28) 



and R abcd = \Rabef£ e ^ cd is the right dual of the Riemann tensor. It will be convenient 
to write explicitly the tetrad components of R* abcd - Since 9 = 0{S), we only need 
R* bcd at 0(9°). Then the non-vanishing components of R* abcd are given by 

2^0123 — -"-0213 — - ft 0312 _ -"1203 ~~ -"1302 _ ~^ -"-2301 — ~3 u \ a )■ 

(11-29) 

Although we do not need them, we note that the following components are not 
identically zero but are of 0(6). 

^1212) ^1313 1 -^1010' -^2323 > -^2020 > an ^ -^3030 ■ 

Further, we may set = in the equations of motion (11-27). 
11.2.1. Lowest Order in S 

We first solve the equations of motion for a circular orbit at r = ro at the lowest 
order in S. For notational simplicity, we omit the suffix of ro in the following. 
For the class of orbits we have assumed, we have v 1 = and v 2 = 0(9). Then the 
non-trivial equations are 



d 

d^ 



1 - Av° + Bv 3 = 0, (11-30) 









H 







|C 2 = 0, ±\ C I I 1 D I [ , l |. (11-31) 



A 


-B 




The equation (11-30) determines the rotation velocity of the orbital motion. By 
setting £ := v 3 /v°, we obtain the equation 

ui + 2w 2 £ + w 3 £ 2 = 0, (11-32) 



Chapter 1 Black Hole Perturbation 



57 



which is solved to give 



iVMr 



A 



(11-33) 



The upper (lower) sign corresponds to the case that v 3 is positive (negative). Then, 
i 1 

0(S 2 ), we find 



with the aid of the normalization condition of the four momentum, v^v^ = — 1 + 



(11-34) 



Note that, in this case, the orbital angular frequency Q v is given by a well known 
formula, 



Qlf r 3 / 2 ± VMa ' 
On the other hand, the equations of spin (11-31) are solved to give 



(11-35) 



/ 






( 


a sin(</> + ci) + (3c2 










cos((p + Cl) 


V 


C 3 y 






— Psm(4> + ci) — ac2 



(11-36) 



where 0_, ci and C2 are constants, and 

A 



o 



VB 2 - A 2 




(11-37) 











( asm(j) \ 


C 2 = -Ci , 






l- 


COS (f> 




I 


C 3 y 




[ —j3 sin <j) J 



The supplementary condition v a ^ a = requires that C2 = 0. The condition CaC° = 1 
implies Cl + Cy = 1- Further since the origin of the time r can be chosen arbitrarily, 
we set ci = 0. Thus, we obtain 



(11-38) 



Here, we should note that ft p ^ fi v in general if a ^ or S ^ (see below). 
11.2.2. First order in S 

Having obtained the leading order solution with respect to S, we now turn to 
the equations of motion up to the linear order in S. We assume that the spin vector 
components are expressed in the same form as were in the leading order but consider 
corrections of O(S) to the coefficients a, j3 and Q p . As we have noted, Eq. ( 11-6| ) tells 
us that v a can be identified with u a to 0(S). In order to write down the equations 
of motion up to the linear order in S, we need the explicit form of R a , which can be 
evaluated by using the knowledge of the lowest order solution. They are given as 



R° = R 3 = 0(6), 
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M 



R 1 = 3 ^rv°v 3 ( 2 + 0(6), 



M 



r 2 = 3jAV + o(e). 



(11-39) 



First we consider the orbital equations of motion. With the assumption that 
v 1 = and v 2 = O(0), the non-trivial equations of the orbital motion are 



v l = Av° + Bv 3 - SR 1 = 0, 
v 2 = (Cv° + Dv 3 )9 - SR 2 . 



(11-40) 
(11-41) 



The first equation gives the rotation velocity as before, while the second equation 
determines the motion in the ^-direction. 

Again using the variable £ = v 3 /v°, Eq. (11-40) is rewritten as 



wi + 2w 2 £ + w 3 e 2 + 3^^£ = 0, 



where S± := S£±. The solution of this equation is 



A 



2r 3/2 



(11-42) 



(11-43) 



Using the relations (11-34), it immediately gives v° and v 3 . From the definition of 
the tetrad, we have the following relations, 



smi 



dt . 9 dip 

asm 2 9 -f- 

dr or 



-a~T + {r +a) — 



(11-44) 



Thus, the orbital angular velocity observed at infinity is calculated to be 



dip 



n<p ' dt r 2 + a 2 + a£VA 



+ O{0 2 ) 



- 3 / 2 ± aVM 



1 



3S ± ±VMr_ - a 
2 r 2 ± aV Mr 



+ 0(9 2 ). (11-45) 



In order to solve the second equation (11-41), we note that v 2 

M 1 + 2£ 2 



Cv° + Dv 3 



r 2 i~e 



+ 0(S). 



Then we find that Eq. (11-41) reduces to 



Ml + 2£ 2 ~ S\\M £ 
^Y^t 2 r3 l-£2 COSl 



EQ ~ r9 and 



(11-46) 



(11-47) 
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where S\\ = SC\\ ■ This equation can be solved easily by setting 9 = 9q cos <f>. Recalling 
that Ql = M/r 3 + Q(S), we obtain 



5|| 



(11-48) 



Thus we see that the orbit will remain in the equatorial plane if Su = 0, but deviates 
from it if Si I ^ 0. We note that there exists a degree of freedom to add a homogeneous 



solution of Eq. (11-47), whose frequency, Qq 



m i • -ie 

l _ £2 



is different from f2 n 



and which corresponds to giving a small inclination angle to the orbit, indifferent 
to the spin. Here, we only consider the case when this homogeneous solution to 9 
is zero, i.e., those orbits which would be on the equatorial plane if the spin were 
zero. Schematically speaking, the orbits under consideration are those with the total 
angular momentum J being parallel to the z-direction, which is sum of the orbital 
and spin angular momentum J = L + S (see Fig. 3)). 




Fig. 3. A schematic picture of the precession of orbit and spin vector, to the leading order in S. 
The vector J represents the total angular momentum of the particle. The vector L is orthogonal 
to the orbital plane and reduces to the orbital angular momentum in the Newtonian limit. In 
the relativistic case, however, these vectors should not be regarded as well-defined. 

Next we consider the evolution of the spin vector. To the linear order in S, the 
equations to be solved are 

(° = AC 1 + C( 2 9 - Sv°C a R a , 
C 1 = A(° + B( 3 + E( 2 , 



( 2 = (C(° + £>C 3 ) 
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( 3 = -BC 1 - D( 2 6 - Sv 3 ( a R a . 
The third equation is written down explicitly as 



c 



i2 



-0(\\nsm(j)cos(j), 



where 



Thus we find 



k := aD — (5C — QpUj^r. 



(11-49) 

(11-50) 
(11-51) 
(11-52) 



Since the spin vector S a is itself of O(S) already, the effect of the second term is 
always unimportant as long as we neglect corrections of 0(S 2 ) to the orbit. 

The remaining three equations determine a, (5 and fi p . Corrections of O(S) 
to a and (3 are less interesting because they remain to be small however long the 
time passes. On the other hand, the correction to ft p will cause a big effect after 
a sufficiently long lapse of time because it appears in the combination of Q p t. The 
small phase correction will be accumulated to become large. Hence, we solve fi p to 
the next leading order. Eliminating £° and £ 3 from these three equations, we obtain 



{B 2 -A 2 )-Q 2 



S±_ / AC -BP 
T V r 



- 



Then after a straightforward calculation, we find 



ni 



M 



1 



3S± ±VM (2r 2 - 3Mr + a 2 ) + ar l l 2 {M - 



-3/2 



- 2 - 3Mr ± 2aVMr 



(11-53) 



(11-54) 



As noted above, f2 p / Q v for S± / 0. The difference Q p — gives the angular 
velocity of the precession of the spin vector, as depicted in Fig. 3. 

11.3. Gravitational waves and energy loss rate 

We now proceed to the calculation of the source terms in the Teukolsky equation 
and evaluate the gravitational wave flux. For this purpose, we must write down the 
expression of the energy momentum tensor of the spinning particle explicitly. We 
rewrite the tetrad components of the energy momentum tensor in the following way: 



-iab 



5^(x-z(t)) 



5^(x-z(t)) 



fi J dr< 



u^+^ a v b ^S dc 



^S b)d v c 



Z(T 




+ 



^=da (B^ b 5 4 {. 



X 



Z\T 




(11-55) 
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The last line gives the definition of A ab and B^ ab . Then the source term of the 
Teukolsky equation is given by Eq. ( 2-14 ) with Eqs. (2-15). 

As we will see shortly, the terms proportional to Su in the energy momentum 
tensor do not contribute to the energy or angular momentum fluxes at linear order 
in S. In other words, the energy and angular momentum fluxes are the same for all 
orbits having the same S±. Thus, we ignore these terms in the following discussion. 
Further we recall that the particle can stay in the equatorial plane if S\\ = 0. Hence 
we fix 8 = tt/2 in the following calculations. 

Using the formula ( |2T3| ), we obtain the amplitude of gravitational waves at 
infinity as 

^Imuj — ^tmuj ' ^Imuj "+" "tmui ^ 11-00 j 

where 

-l 

A r . 



J £muj 



(.mui 







i^Kn + irnB% n - B r nn — 



-t. 



L 2 p^ 2 S] 



t„3 



Cm 



, . ryin 
u±S £rmj 



5 {uj — mf2) 



6=n/2 
dt^- 1 r 

dr 



1 

rA 

A _ 



L £mui 



r=ro 
t 



d 



J | QOW 



e=n/2 \[A 



' d 2iK 4 
dr A r 



L £muj 



-1 r 



A,. 



iuiBjf.jf, + imB mm B, 



d_ 

"' m dr 



x (-2-5^)0=^/2 



dr 2 



1 %K\ d (%K 
2 ! " + ~A J dr~ ~ \ A~ 



+ 



2iK K 2 



Ar A 2 



L £muj 

r=r 

(11-57) 



and 



B^- 



B» - 



1 1 



41 ~e 



{l-S ± ((2u 1 +u 3 )Z + u 2 )} 
1 / r 2 + a 2 t rr " 



4r 1-e 



e + a , -VZf, 0, -w + i 



r 2 + a 2 



1 1 



h Sj -i^ (°' ^ °' °) • 



(11-58) 
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The Lorentz factor dt/dr which appears in Eqs. (11-57) is calculated from Eqs. (11-44) 

as 

dt 1 / r 2 + a 2 x 



a£ + — j=-\. (H-59) 



In general, as we have seen in the preceding sections, when the orbit is quasi- 
periodic the Fourier components of gravitational waves will have a discrete spectrum; 



Zgmuj = Vi(w- u n )Ze muJn . (11-60) 



n 



Then the time-averaged energy flux and the z-component of the angular momentum 
flux are given by the formulas ( 2-3l| ) and ( 2-32| ), respectively. In the present case, 



since we may regard the orbits to be on the equatorial plane, the index n degenerates 
to the angular index m and u n is simply given by mfl^ (n = m). Hence we eliminate 
the index n in the following discussion. Here we mention the effect of nonzero S\\ . If 
we recall that all the terms which are proportional to Su have the time dependence 
of e ±J ^p T , we find that they give the contribution to the side bands. That is to say, 
their contributions in Zi muJ are all proportional to S(u — mfi ± Op). Then, since the 
energy and angular momentum fluxes are quadratic in Zi mu}n , they are not affected 
by the presence of Sn as long as we are working only up to linear order in S. 

As before, in order to express the post-Newtonian corrections to the energy flux, 
we define 7j£ muJ as 

'dE\ 1 fdE s 



u , o , „ , Vim, (11-61) 
at J lm 1 V at / N 



where (dE/dt)^ is the Newtonian quadrupole formula defined by Eq. ( 4-19[ ). 



We calculate r]£ m up to 2.5PN order. Keeping the 5-dependent terms, the results 

are 

W - ( 1 ».3 1 „„,4 535 .,5 



^=\-T2 V -8 qV " 1008^ S ' 
W _ 10935 5 . 



^±3 



896 



V s, 



(s) 20 5 „ 
r?Qj_o = — v s, 

/3±2 g 3 , 

= 7^ 5 -s, (11-62) 

/3±1 8064 ' v ; 

where q = a/M and s := S±/M. The rest of rjitl are all of higher order. We should 
mention that if we regard the spinning particle as a model of a black hole or neutron 
star, S is of order fi. Therefore the correction due to S is small compared with the 
/5-independent terms in the test particle limit n/M <C 1. 
Putting all together, we obtain 



dE\ _ /dE\ 
~dt I ~ Vdt) 



N 



„ 1247 2 / 73 25 A , 

1 v z + 4-7T q s v 6 

336 V 12 4 



Chapter 1 Black Hole Perturbation 63 



44711 33 



71 

+ —q's I v 



+ 



9072 + 16' 



8191 3749 2403 

-7T + ^TT? + -TT^S 



672 



336 



112 



(11-63) 



Since v is defined in terms of the coordinate radius of the orbit, the expansion with 
respect to v does not have a clear gauge-invariant meaning. In particular, for the 
purpose of the comparison with the standard post-Newtonian calculations it is better 
to write the result by means of the angular velocity observed at infinity. Using the 
post-Newtonian expansion of Eq. (11-45) 



1 



v 3 + -qsv + 



o (»«)), 



(11-64) 



Eq. (11-63) can be rewritten as 



dE\ 
~dt) 



dE\ 
~dT)^ 



1247 



+ 



+ 



336 
44711 



X 2 + 47T 



11 £ 
T 9 ~ 4 S 1 ''" 



33 o 31 ^ \ 
9072" + 16 9 +Y qs]x 



8191 59 

7T -| ( 

672 16 



13 

16' 



s x 
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where x = (MJ?^) 1 / 3 and (dE/dt) N is the Newtonian quadrupole formula expressed 
in terms of x, Eq. ( |6-6[ ). Since there is no sideband contribution in the present case, 
the angular momentum flux is simply given by (dJ z /dt) = Q~ l (dE / 'dt) gw • The 
result (ll-65l_is consistent with the one obtained by the standard post-Newtonian 
approach eJ)'ES to the 2PN order in the limit fi/M — ► 0. The s-dependent term of 
ojder x 5 is the one which is newly obtained by the black hole perturbation approach 



§12. Black hole absorption 



When a particle moves around a Kerr black hole, it radiates gravitational waves. 
Some of those waves are absorbed by the black hole. We calculate such absorption of 
gravitational waves induced by a particle of mass fi in circular orbit on the equatorial 
plane around a Kerr black hole of Mass M. 

The post-Newtonian approximation of the absorption of gravitational waves into 
the black hole horizon was first calculated by Poisson and Sasaki in the case when a 
test particle is in a circular orbit around a Schwarzschild black holeBiP. In this case, 
the effect of the black hole absorption is found to appear at 0(v 8 ) compared to the 
flux emitted to infinity and it turns out to be negligible for the orbital evolution of 
coalescing compact binaries in the near future laser interferometer's band. On the 
other hand, the black hole absorption appears at 0(v 5 ) if a black hole is rotating. 
That calculation was done by Tagoshi, Mano, and TakasugiES. 
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In order to calculate the post-Newtonian expansion jjf ingoing gravitational 
waves into a Schwarzschild black hole, Poisson and SasakiKIP used two type of rep- 
resentations of a solution of the homogeneous Teukolsky equation. One is expressed 
in terms of the spherical Bessel functions which can be used at large radius, and 
the other is expressed in terms of a hypergeometric function which can be used near 
the horizon. Then two type of expressions are matched at some region where both 
formulas can be applied. They got formulas for a solution of the Teukolsky equation 
which can be used to calculate ingoing gravitational waves to 0(i> 13 ), although they 
gave formulas for ingoing waves only to 0(v 8 ). 

Here, we first review the method found by Mano, Suzuki and Takasugi , since 
it is the only existing method by which higher order post-Newtonian terms of the 
gravitational waves absorbed into a rotating black hole can be calculated. We note 
that this method is also the only existing method that can be used to calculate the 
gravitational waves emitted to infinity to an arbitrarily high post-Newtonian order. 
Then we calculate the energy flux absorbed into the horizon to 0(u 13 ), i.e., 0(v 8 ) 
beyond the lowest order flux absorbed into the horizon, for circular orbits on the 
equatorial plane of a Kerr black hole. 

12.1. Analytic solutions of the homogeneous Teukolsky equation 

Analytic series solutionsof the homogeneous Teukolsky equation were found by 
Mano, Suzuki and Takasjugic 3 ! and various properties of the solution were discussed 
by Mano and TakasugiS. Here, we follow the notation of Ref. 55) except that 
we focus on the case of spin weight s = —2. In this method, the solution of the 
radial Teukolsky equation (2-3) is represented by two kinds of expansion. One is 



given by a series of hypergeometric functions and the other by a series of Coulomb 
wave functions. The former is convergent at horizon and the latter at infinity. Then 
the matching of these two solutions are done exactly in the overlapping region of 
convergence. 

First we consider the solution expressed in terms of hypergeometric functions. 
The solution which satisfies the ingoing wave boundary condition at horizon is ex- 
pressed as 



oo 

X 

n=— oo 



2_. a^F(n + v + 1 — ir,—n — u — ir,3 — 2ie + ;x), (12T) 



where F(a, b, c, x) is the hypergeometric function and 

w(r " r+) - TT^, T =t^, e± = l±I. (12-2) 



en v k 

The coefficients a v n obey a three terms recurrence relation, 

«X + i + « + 7X-i = 0> (12-3) 

where 

v ien(n + v — 1 + ie){n + v — 1 — ie){n + v + 1 + it) 



{n + v+ l)(2n + 2i/ + 3) 
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(in = -A - 2 + (n + v){n + v + 1) + e 2 + e(e - mg) 
e(e - mg)(4 + e 2 ) 



+ 



{n + v){n + v + I) 



v _ ieK(n + v + 2 + ie)(n + v + 2 — ie)(n + v - ir) , . 

7 " " (n + ^)(2n + 2z^- 1) ' 1 ' j 

The series converges if z/ satisfies the equation, 

Rn^)L n ^) = 1, (12-5) 
where R n {v) and L n (v) are the continued fractions defined by 

jf fjA — u n In 

a v a u 

""n 

The range of convergence is < (— x) < oo for physical x. We can prove that if v is 
a solution to Eq. (12-5), then so is —v — 1. Since we can set n in Eq. (12-5) to an 
arbitrary integer, a convenient choice is to put n = 1. Further, for convenience, we 
may set a$ = Cq = 1. It then follows that we have 

aZ v - l = a v n . (12-7) 

This implies R inu = ff 11 ^ 1 = R in . Consequently, for (-x) > 0, Eq. (12-1) can be 
rewritten as 

R m = e ieK {Rl + R^ l ) 1 (12-8) 



a n+l Pn + ln^n-\\V) 



where 



r(3-2ie+)r(2n + 2i^ + l) 



oo 

X 

n=— oo 



E 



r{n + u + \ — iT)r(n + v + 2> — ie) 



where 



xx n F(-n -v-iT,-n-v + 2-ie, -2n - 2u; 4), (12-9) 



uj(r — r) 

x=— '-. (12-10) 

en 



Since v — {—v — 1) = 2v + 1 is not an integer in general, the solutions Rq and Rq u1 
form a pair of independent solutions. 

The other series solution which is convergent at infinity is expressed in terms of 
the Coulomb wave functions;^ 

R o - 4 - t f " + Jj-*'" j^f^<^>. (mi) 

where z = u>(r — r_) = e/tx, (a) n = -T(n + a)/r{a), and F n+U (z) is the Coulomb 
wave function given by 

F n+u{z ) = e-{2 Z )^ z r{n + V + * + K) 



r(2n + 2v + 2) 



GG 



x$(n + v + 3 + ie, 2n + 2v + 2; 2iz), 

where &(a,b,z) is the regular confluent hypergeometric function 0*. A crucial ob- 
servation made by Mano, Suzuki and TakasugiH is that the coefficients a y n obey the 
same recurrence relation as that for the hypergeometric type solution, Eq. ( 12-3| ). 
The series (12-11) converges in the range z > en if v is a solution of Eq. ( |12-5| ). The 
solution Rq can be decomposed into a pair of solutions, a purely incoming wave at 
infinity RK and a purely outgoing wave at infinity R v _. Explicitly, we have 



where 



R V C = R V + + R U _, (12-12) 



R-- 2 » e -^(»+z) F ^ ± 3 ± ie l e -i^ u+ie + (z - e K ) 2 " ie + 
F[y — 1 — ie) 



x J2 i n an(2z) n y(n + v + 3 + ie,2n + 2is + 2;2iz), (12-13) 

n=— oo 

oo 

R u_ =2 u e -ne e -m(u-l) e i~z~u+ie + ^_ e ^2-ie + g -n 

n=— oo 

<(2z) n ^(re + v - 1 - ie, 2n + 2^ + 2; -2iz), (12-14) 



(i/ - 1 - ie) n v 



(v + 3 + ie) n 

where ^(a, b, z) is the irregular confluent hypergeometric function©. By definition, 
the upgoing solution R up is given by 

R u p = R u _. (12-15) 

We see that the above two kinds of solutions are convergent in the common region 
1 < x < 00. Then comparing the asymptotic behaviors of Rq and R v c for x — > 00, 
we find that they have the same characteristic exponent, ~ x v , hence describe the 
same solution up to the normalization factor. Therefore by comparing each power 
of x we have 

R V Q =K V R V C , (12-16) 



where 



(2eK)-^ 2 - r 2 2 iT(3 - 2ie+)r(r + 2u+ l)r(r + 2v + 2) 

Tiv + 1 - ir)T(v + 3 - ie)F{f + v + 3 + ie) 
r{v + 1 + ir)r(u - 1 + ie) 



r{r + v + 1 + fr)r(f + v - 1 + ie) 
(™ (f + 2v + l) n (u-l-ie) n v 
X L% (n-f)! (v + 3 + *e) n °* 
/ ? 



x I V , J=£ -|- 1 - i t< ■ (12-17) 

where f can be any integer and K v is independent of the choice of f. 
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The gravitational wave absorbed into the black hole is expressed by Eq. Q2T2|) . 
Hence we need to know the amplitudes B inc and B tians of R in and C trans of R up , 
defined in Eq. (|2-7j ). The asymptotic ingoing amplitude at horizon, £? trans ; Q f R m is 
readily obtained from Eq. (12-1) as 



B l 



CK 



2,s 



In 



K E < 



(12-18) 



Similarly, the asymptotic outgoing amplitude at infinity, C ,trans , of R up is obtained 
from Eq. (12-14) as 



£<trans 2l+* e g — 7re /2g— l) [i ,3gielne 

E {u — 1 — ie) n 
,_ „ (V + 3 + IfL 



(12-19) 



On the other hand, a bit of work is necessary to obtain the asymptotic incoming 
amplitude at infinity, B mc , of R m . Setting the asymptotic behavior of Wq at z — > oo 

as 

R u c -> A u + z- 1 e-^ z+elnz ^> + A u _z 3 e^ z+€lnz \ (12-20) 

we find 

A „ = 2 -3-i £e i(7r/2)(i/+3) e -7re/2£> + l+if) y 

r(u — 1 — ie) ^ n ' 

v ' n=— oo 

^ =2 i+- e -^/2)(^i) e -W2 ^ (~l) n ^~^~^" <. (12-21) 

(z^ "I - 3 ~\~ %e)rt 

n=— oo v 

Because of Eq. (|l2~^ ), A u ± and are related to each other as 



-v-l 



,smn(u + ie) 



-ie 



sin7r(i/ — ze) 



4^ 



A!"' 1 = ie inu A u _ . 



(12-22) 



With the help of the above relations, we find from Eqs. (12-8) and ( |12-16| ) the 
asymptotic amplitudes of R m at infinity as 



B 1 



K v — ie 



, shimf + ie) 



sm-7r \v — ie) 



v-l 



v —ie In e 



Ale 



B ref = jenj, 



K v + ie nv K. 



v-l 



v ie In f 



A v _e 



(12-23) 



So far our discussion has been on exact analytic series expressions for the homo- 
geneous Teukolsky functions. Now we consider their post-Minkowski expansion by 
assuming e«l. Provided we set Oq = 1, we see from Eqs. (12-4) that a^, 7^ = 0(e) 
and (3^ = 0(1) unless the value of v is such that the denominator in the expression 
of or 7^ happens to vanish or (3% happens to vanish in the limit e — > 0. Except for 
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such an exceptional case, it is easy to see from Eq. ( |12-3| ) that the order of a v n in e in- 
creases as \n\ increases. Thus the series solution naturally gives the post-Minkowski 
expansion. 

For the moment, let us assume that the above mentioned exceptional case does 
not happen for n = 0. Then we have R\[y) = 0(e) and Lq(u) = 0(l/e). This 
implies (3q + JqL-i(v) = 0(e 2 ). Then assuming L-\(y) = 0(e), we must have 
(3q = 0(e 2 ). Using the expansion of A given by Eq. (|3-1|), we then find v = £ + 0(e 2 ) 
or v = — i — 1 + 0(e 2 ). Since we know that —v — 1 is a solution if v is so, we may 
take the solution v = I + 0(e 2 ) without loss of generality. Then the assumptions 
that R\(v) = 1/Lq(v) = 0(e) and L_i(i/) = 0(e) are justified. Further it is easily 
seen that Rn(v) = 0(e) for all n > 0. On the other hand, for n < 0, a u n = 0(1) at 



n = - 



x\y. 



1 and 0% = 0(e 2 ) at n = -21 - 1. 
= 0(1/ e). To summarize, we have 



Thus we have L_^_i(V) = 0(1) and 



R n (u) 



0(e) for all n > 0, 



-i 



0(1) 



-H-l 



-2( 



LAv) 



0(e) for all the other n < 0. 



O(IA), 



(12-24) 



With the above results, the post-Minkowski expansion of the homogeneous Teukolsky 



functions can be obtained with arbitrary accuracy by solving Eq. ( 12-5 ) to a desired 
order and by summing up the terms to a sufficiently large |n|. For our present 
purpose, we need v which is accurate to 0(e 2 ). Solving Eq. ( 12-5 ) to this order, we 
find 



+ 



21 + 1 



+ 1 



+ 



_ (l-2) 2 (l + 2) 2 
(21 + \)(2l + 2)(2t + 3) ~ (2i - l)2£(2£ + 1) 



(£-l) 2 (i + Zf 



(12-25) 



Interestingly, v is found to be independent of the azimuthal eigenvalue m to 0(e 2 ). 

The post-Newtonian expansion in the near zone is given by further assuming 
e -C z <C 1 in the series solution ( 1 2 - 1 1| ) and expand it in powers z. For evaluation 
of the black hole absorption, we need the post-Newtonian expansion of R up which 
is obtained from Eq. (12-14). The explicit post-Newtonian formula for i? up and the 
asymptotic amplitudes B mc , £? trans and (7 trans to 0(e 2 ) are given in Appendix [H|. 

12.2. Absorption rate to 0(v 8 ) 

In this subsection, we evaluate the energy absorption rate by a black hole. The 



energy flux formula is given by Teukolsky and PressEf as 



fdE h 
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V dtdQ 
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duJ 



2S 2 



(111 



128ujk(k 2 + 4e 2 )(A; 2 + 16e 2 )(2Mr_ 



2tt 



\C\' 



1 7 R I 2 



(12-26) 
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where e = «/(4r+) and 



A + 2) 2 + 4awm - 4a 2 w 2 J ^A 2 + 36aum - 36a 2 w 2 
+(2A + 3)(96a 2 w 2 - 48aum) + 144w 2 (M 2 - a 2 ). 



(12-27) 



The calculation of Z^ mijJ is parallel to the calculation of Z^^ except that R m 
is replaced by R up . The solution of the geodesic equations are given in section ||. 
Using that solution, we have the amplitude of the Teukolsky function at the horizon 
in Eq. ([M2|) as 



J lmu) 



2TrB tiaBB S(u - mf2) 
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? up 



^£muj{ J ^ nn -^-rnnO ^mmo} 



■ l-^m til T ml; T 



2 D U P 



d 2 R 



dr 2 



mm2 



r=ro,d=TT /2 



dr 

= «J(£j - mfl)Zf m . (12-28) 
From Eq. ( p^ffi ) and (12-28), the time averaged energy absorption rate becomes 
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uj=mQ 



(12-29) 



As in the case of the Teukolsky function at infinity, we can show that Z?_ 

= {-lYzf muJ . Then, from Eq. (12-29), we have (dE/dt)i t - m = (dE/dt)i, m . 

In order to express the post-Newtonian corrections to the black hole absorption, 
we define rjf m as 



dE 
~dt 



dt J 



(12-30) 



N 



where (dE/dt)^ is the Newtonian quadrupole luminosity at infinity, Eq. ( 4-19| ). In 
Appendix |, we show rjg m . 

The total absorption rate to 0(v 8 ) beyond the lowest order is given by 
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where 



A r . 



^ (0, ( 3+ tS?W 0, ( 3 



rag 



Chapter 1 Black Hole Perturbation 



71 



B n 



2i 



2 



V> (0) (^3 + 



mq 



(0) 



0~ 



rag 



+ 3 



mg 



mg 



(12-32) 



and ip( n > (z) is the polygamma function. We see that the absorption effect starts at 
0(v 5 ) beyond the quadrapole formula in the case q ^ 0, while for q = 0, the above 
formula reduced to 

DHfUw). < 12 -> 

as was found by Poisson and Sasaki©. We note that the leading terms in (dE/dt)n 
are negative for q > 0, i.e., the black hole loses the energy if the particle is corotating. 
This is because of the superradiance for modes with k < 0. In Appendix [|, we also 
show (dE/dt)n written in terms of x = (Mfitp) 1 ' 3 . 

It is not manifest from Eq. (12-31) that if it has a finite limit for \q\ — > 1. But 
by using the formulas, 
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lim ip {k) 3 + 
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we obtain the limit of {dE/dt)n as 



0, (k + 0), 



(12-35) 



<A±i ( dt ) H 



N 



9 49 « 2 , , , 

- — —V + 4 7T + 



16 Q 



133 
T2" 



6817 q 4 /535 

TT ^ + 

336 9 V 16 



97 



Z3424 g , 1712 g 3424 g , , . 

+ I ttzt A ln(2) + —— 7 + — — ^ ln(u) 



V 105 g| 



105 



g| 105 |g 



+ 



3647533 q 289 g 2 g , , 

— — 7T — 16/3 7T — \ V 

22050 |g| 6 |g| ' 
/ 84955 55873 



V 336 



+ 



672 



7T V 



( 14077 g , 16441 g 34987 g , , , 

V 60 |g| 140 |g| 210 |g| 



4057965601 g 1289 g > s 

7T ] V 



193 2 _q^ 
12 *" |g| 



6350400 



9 |g| 



(12-36) 



Appendix A 

Spheroidal harmonics 



72 



In this Appendix, we describe the expansion of the spheroidal harmonics -2<S^ 
to 0{{auj) 2 ). 

The spheroidal harmonics of spin weight s = — 2 obey the equation, 

1 d ( n d) 99 9 „ (m — 2cos#) 2 
\ sin 0— ^ - orur sm 2 6> - 



sin 6» d6» I d(9 J sin 2 9 

+4aw cos - 2 + 2maw + A _ 2 5g£ = 0. (A-l) 

We expand - 2 S^ and A as 

_ 2 Sg£ = _ 2 P< ro + au;S« + (au,) 2 ^ + 0((a W ) 3 ), 

A = Ao + awAi + a 2 uj 2 \2 + 0((aw) 3 ), (A-2) 

where - 2 Pe m are the spherical harmonics of spin weight s = —2. We set the normal- 
izations of -iPlm an d -2'S'^m as 

I* \-2Pe m \ 2 sm6d6 = F |_ 2 ^| 2 sinfldfl = 1. (A-3) 

JO JO 

Inserting Eq. (A-2) into Eq. (A-l) and collecting the terms of the same order to 
(aw) 2 , we obtain 

[Co + A ] - 2 Pim = , (A-4) 
[Co + A ] Sj& = -(icosO + 2m + Ai) _ 2 ^ m (A-5) 
[Co + A ] 5^ = -(4 cos 6 + 2m + Aijsgj - (A 2 - sin 2 6) - 2 P lm , (A-6) 

where 

Id/ n d\ (m-2cos6>) 2 , A „, 

sin 6 dd \ dd J sir 6* 

The lowest order equation (A-4) says we have Ao = {£ — 1)(£ + 2). 

The first order correction to the eigenvalue, Ai, is obtained by multiplying 
Eq. (A-5) by - 2 Pg m from the left hand side and integrating it over 0. The result is 



To obtain S^, we set 



Sffi = E -2i>m - (A-9) 



We insert this into Eq. (A-5), multiply it by - 2 Pe m and integrate it over 6. Then 
noting the normalization of the spheroidal harmonics, we have 



4m=\ {£' -!){£' + 2) -{£-!){£ + 2) 

o, £' = e. 



J d(cOS 9) cos 6 - 2 Pe m , £'¥=?, 
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Hence cf is non-zero only for £' = I ± 1, and we obtain 



J+l 

4m 



+ 3)(£- l)(£ + m + l)(£-m + l)] 1 / 2 



(£ + l) 2 [ (2^ + l)(2£ + 3) 

(£ + 2)(£ -2)(£ + m){£ - m)l 1/2 



(2* + l)(2£- 1) 

The next order equation can be solved similarly. The second order correction 
to the eigenvalue, A2, is obtained by multiplying Eq. (A-6) by - 2 Pt m from the left 
hand side and integrating it over 9. We find 



A 2 = -4 J d(cos &)-%Pi m cos 9 sf^ + J d(cos 9)^ 2 P hn sin 2 9 _ 2 P fe 

= -2(1 + 1X4+ 1 ) 2 + 2£(c i hn 1 ) 2 + 1 - J d(cos 0)-2^m cos 2 9 - 2 PdA-W) 
where the last integral becomes 
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As before, to obtain 5) , we set 



sfi = E4» - 2 ^ 



(A-ll) 



Inserting Eqs. (A-9) and (A-ll) into Eq. (A-6), multiplying it by -. 2 Pl'm an d integrate 
it over 9, we obtain 
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where < ji, j 2 , mi, m 2 | J, M > is a Clebsch-Gordan coefficient. For 
non- vanishing (£' / £) are given explicitly as 



2 and 3, the 
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4+ 2 = ^^(3 - m) 1 /2 (3 + m) V2 (4 _ m) i/ 2 ( 4 + m ) 1/2 , 



for £ = 2, and 
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for I = 3. As for d^ m , it is determined by the normalization of _2<S'|^, i.e., 
1 = J d(cos6)\- 2 S£m\ 2 
= / d(cos0) <^ (- 2 P| m ) 2 + 2aw^ 

+ ( a ^) 2 X^m c £m-2-P£'m-2-P<?"m 
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Then we have 
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The operators , and 



(A-13) 



In this Appendix, we show the operators for n = 2, 3 and 4 which appear 
in Eq. (ggof). 
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Appendix C 

4PN formulas for R t 




in 

tm 

In this Appendix, we show the post-Newtonian expansion of R m in the near zone, 
where z = u>r <C 1, for a Kerr black hole which are needed to evaluate gravitational 
waves at infinity to 0(v 8 ). For convenience, we recover the indices im on R m and 
give the formulas for ujR 1 ^. 
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The ingoing Reg ge- Wheeler functions to 0(e' 



In this Appendix, we present a method to calculate the ingoing Regge- Wheeler 
functions to 0(e 3 ) which are needed to calculate the luminosity to 0(v n ) beyond 
Newtonian in the Schwarzschild case. In the Schwarzschild limit, q — > 0, and solving 



Eq. ( 3-25 ) recursively is in principle straightforward. Since the general homogeneous 
solution to the left-hand side of it is given by a linear combination of the spherical 
Bessel functions je and ni, one can immediately write the integral expression for 
Q n ^ . Noting that jgn\ 



nd'i 



(n) 



1/z 2 , we have 



where the source term in the Schwarzschild case is given by 

1 d 



W (n) = ^(1)^-1)] 



,d_ 

dz 



z 3 dz 



(D-l) 



(D-2) 



We perform the above indefinite integral and set the appropriate boundary condition 
by examining the asymptotic behavior at z — > order by order. 



D.l. General remarks 

(n) 

Let us first consider the boundary conditions of Q . In the Schwarzschild case, 
the original Regge- Wheeler ingoing wave function Xf 1 is related to ^f 1 as 



X 



in —ie\n(z—e)An 

■ ze q> 



(D-3) 



and it has the asymptotic behavior given by Eq. ( 3-12 ), which in the present case 
reduces to 



XI 



(Af 
U 



+ Af c e~ 



trans _,— iojr* 



for Z* 
for Z* 



oo 



-oo. 



(D-4) 



Thus, noting that z* = z + eln(z — e), the boundary condition of ^f 1 is that it is 
regular at z* — > —oo (z — > e). To implement this boundary condition to £} n , we need 
a different series expansion of it; a series in terms of the variable x := r/2M = z/e 
around x = 1. We write this expansion as 



d 0} (x) + 6^ i| (x) + 6 2 er i (^) + 



;{!} 



2 C {2}, 



(D-5) 



On the other hand, there are two independent solutions expanded by functions writ- 
ten in terms of z. We denote the solution whose zeroth order is given by ji (ni) as £j£ 
(ini)- Then the general solution is given by Cj£ji + c n ^ n g. As can be shown by using 
the result of Poisson and Sasaki©, ^j ^ does not have terms matched with £ n £. A 
term matched with £ n £ first appears from ^j 1 ^. For x — > oo (<-> e <C z), cj°\ x ) De_ 
haves as x e = z e e~ e ~ e~~ On the other hand, e^j 1 ^ contains the term that behaves 
as ex~ i_1 = e i + 2 z ~ i ~ 1 ~ ^ +2 Cni- Therefore, in the sense of the post-Minkowskian 
expansion, the inner boundary condition affects the ingoing wave solution at and 
beyond 0(e 2i+2 ). However, in the post-Newtonian sense, since we evaluate ^ in the 



so 



near zone, i.e., for z = 0(v), the contribution from becomes 0(v u+s ) relative to 
Since t > 2, we find that the inner boundarvxondition affects the homogeneous 

solution at and beyond 0(v 13 ) in the near zone.tH 

Since j t = 0{z e ) as z -»■ 0, we have Xf -► C^+^e"***, or ^ rans = 0(e e+1 ). 

On the other hand, from the asymptotic behavior of je at z = oo, we find A^ 110 and 

^ ef are of order unity. Then using the Wronskian argument, we obtain 

I A trans 1 2 

\Af*\ ~ \Af\ = , -J 7 ,' f , = 0(e 2E + 2 ). (D-6) 

|y^, nc | + \Af\ ^ J V 7 

Thus | A^ nc | = |^ ef | until we go to 0(e 2i+2 ) or more. This fact implies that we can 
make A^ c and A r ^ to be complex conjugate to each other to 0(e 2i+1 ). Hence the 
imaginary part of Xf 1 , which reflects the boundary condition at horizon, appears 
at 0(e 2l+2 ) because the Regge- Wheeler equation is real. This is consistent with the 
argument given in the above paragraph. Provided we choose the phase of Xf 1 in this 

way, Im J for a given n < 2£ + 1 is completely determined in terms of Re 



for r < n — 1. 

is explicitly, let us < 

; (n) _ An) . (n) 



(n) 

To see this explicitly, let us decompose Q into the real and imaginary parts: 



Inserting this expression into Eq. ( D-3 ) and expanding the result with respect to e 
by assuming z 3> e, we find 

XT = e-^z (j t + e(/W + + ^(/f + fcf>) + £ 3 (/ f) + ig f)) + . . .) 
* | + + e 2 (/f + ff « In z - ^(ln z) : 



+« f^P " j/ In z) + e 2 (pf + -Jt - In 
+e^gf -\ (In z) 2 gV- In zfP + i/« 



Hence we must have 



7| — Jt in 2 , <fy — - -j£ + J e mz, 
,(3)_ f I(i-^3 JL^. , 1.^(2) 



lOn^-^^-i/^+lnx/f, (D-9) 



*' In Ref. [Tl|), it was erroneously argued that the outgoing gravitational waves are unaffected 
by the inner boundary condition until we reach 0(e 6 ) = 0(v is ). As shown here, this is true only in 
the post-Minkowskian sense. 
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(n) 

For completeness, we also give the relation between the functions /) and the con- 
ventional post-Newtonian expansion of Xf 1 : 



x 



ra=0 



X, 



(0) 



X. 



(3) 



zff ] = zh , X? = ^ , XP = z [ff > + ii,(ln z) 
*(/f + 5/i 1) dn,) 2 -^ 



(D-10) 



Now we turn to the asymptotic behavior at z = oo. Let the asymptotic form of 



/i n) be 



fl n) ^ PPji + QT >n t asz^oo (ra = 1,2,3). (D-ll) 

Then noting Eq. (D-9) and the equality e ~ leln ( z - e ) = e~ lz * e lz , the asymptotic form 
of Xf 1 is expressed as 



XI 



\e--*(zhfe-) l + e {p«+*(gf +ln.)} 

V{(pj 2) -gf in.) +*(gf +pU\nz)} 

{ (/f " Of ^ *) + < (Of + Pf In * + i (In ,) 3 ) } + 

{ 



l + e\pW-i(QW-]nz)} 



+e*{(pP+QV\nz)-i(Qf - P^lnz)} 

+e 3 | (pf } + Qf In z) - i (gf - pf } In z - ~(ln zf^j X + -(-D12) 

Using the asymptotic behavior of tip and faf given in Eq. (Hll), the incident 
amplitude ^4™ c can be readily extracted out as 

A nC = V +1 e-- ln£ [l + e {p, (1) + i (gf + In z) } 

{ (^f " Of In *) + i (gf + In *) } 
+e 3 { (pf - Qf In *) + i (gf + Pf } In * + I (In ,) 3 ) } + • • (t-13) 



(D-14) 



where note that the definition of r* , Eq. (2-9), in the limit q — > is 

r - 2M S 



ojr* = u [r + 2M In • 



2M 



— e In e, 



which gives rise to the phase — ie In e of ^U 1 
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An important point to be noted in the above expression for A" 1C is that it contains 
In z-dependent terms. Since A^ c should be constant, and should contain 
appropriate In z-dependent terms which exactly cancel the In z-dependent terms in 
the formula (D-13). 

D.2. Basic formalism for iteration 

Here we derive the formulas necessary to perform the iteration scheme. 
D.2.1. Definitions 

We introduce the following functions, 



B jf= 


[ zjojodz = 
Jo 


2 C ' 




B n j'- = 


/ zn j dz = 
Jo 




Bj n '- = 


/ zj n dz = 
Jo 






B nn '.= 


/ znQnodz = 


— Bjj + lnz, 


(D-15) 



where 

f 2z , s'my f°° , sin-y ir 

S = / dy^- = - \ dy^ + -, 

Jo y Jx y 2 

C= [ 2Z dy C ° Sy - l =- r^- 7 -ln2,, (D-16) 
Jo y ^ y 

and the lower bound z* of the integral for the definition of B nn is adjusted so as to 
make B nn equal to the last expression of the line. 

As an extension of these integral sinusoidal functions, we further introduce the 
following functions: 

Bjj:= j Z dzzjoDft, (D-17) 
B nJ := ["dzznoDg, (D-18) 

J z, 

where J stands for a sequence of j and n, say, J = jnnj, and we have also introduced 
an extension of the spherical Bessel functions by 

Dj:=j e , D?:=n t , (D-19) 

and 

Df := neBjj-jiBnj, 

D\ J := j e Bjj + n e B nJ . (D-20) 

We adopt the following rule to determine the lower bound of the integrals in Eq. (D-18). 
Whenever we can put = 0, we do so, which is always possible when the sequence 
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J ends with j. On the other hand, in the case when J ends with n, there may 
appear in the integrand the square of no which causes logarithmic divergence if we 
set = 0. In such cases, we use the relation, 

4 = ^-ji, (D-21) 

to replace Uq with the right hand side and extract out the logarithmically divergent 
term due to 1/z 2 . Then we set = for the jg term, while we set z* = 1 for the 
1/z 2 term so as to make the resulting logarithmic term zero at z = 1. For J of two 
indices, this is how we have defined B nn in Eq. (D-15). For J of three indices, this 
applies to B n j n . Specifically it is given by 

Bnjn = B j3j - Bjj In z + ~(ln z) 2 . (D-22) 

For convenience, in what follows we call Bj the generalized integral sinusoidal func- 
tions and Df the generalized spherical Bessel functions. 

Note that all the Bj whose J end with n can be expressed in terms of those 
whose J end with j. For example, for J of three indices, we have 

^jj n — -B nn j -\- In zB n j , 

Bjnn — ^Bjjj -^nnj hi ^Bjj , 



Bnnn — ^^njj ^jnj \XLzB n j , (D'23) 



together with Eq. (D-22). Using these relations, we can express all the Df whose J 



end with n in terms of those whose J end with j. 

Further we introduce the following indefinite integral operator, 

F k/ [X] := n k J dzj e X - Jk J dzn e X, (D-24) 

for a function X. Note that Eq. ( p-l| ) is expressed in terms of this operator as 

CP =F e>e [wW]. (D-25) 
We also introduce the following operator, 

H J k [Y] := n k J dzD J Y - j k J dzD^Y, (D-26) 

where Y stands for a linear combination of the generalized Bessel functions with the 
coefficients given by linear combinations of z m {hiz) n (m < 1, n > 0) and Y denotes 
the quantity which is obtained by replacing j, n, D^ J and D nJ with n, j, D nJ and 
Di J , respectively, in the expression of Y. By definition we see that 

F efi [z m (\nz) n j ] = Hl[z m (lnz) n j }. (D-27) 
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D.2.2. Basic formulas 

The spherical Bessel functions satisfy the recursion relation, 

2m + 1 

Cm-1 + Cm+1 = Cm > (D-28) 

where Cm = j m or n m . Note that 

n m = (-ir +1 i-m-i , ; m = (-l) m n_ m _ x . (D-29) 
The same recursion relation holds for the generalizes spherical Bessel functions, 

D^ + D% +1 = 2 -^Di J , (D-30) 

where ( = j or n. Further the relations the same as Eqs. flDggD hold for D^ J , 

D nj = ( _ r+ i jDi J m _ ij D iJ = ( i -)m J ynJ n _ v (D-31) 

The derivative recursion relation for the spherical Bessel functions is 

^6 = i*U-i -{£+ i)0+i} , (d-32) 

and this extends to the generalized spherical Bessel functions as 

-^F = _L_ _ {e + 1)2 ^j _ ^i^. (D . 33) 

Useful integral formulas for the spherical Bessel functions are 

/dz CmCn = 7 w I rTT (CmCn+1 + Cm— lCra) CmCn' 

(m — n)(m + n + 1) m — n 

(m ^ n, — n — 1), 
/ ^^=2^{/^CoCo-^(CoCo+2ECmC + C^|, 

dz zCmCn = dz ^Cm-lCn-l j (Cm— lCn-1 + CmCn) ; 

J m + n 

dz zQCi = f dz <oCo - y { Co Co + E Q + CmC + ](iQ \ , 

I m=l ^ ' J 

/" zQCl-i = J dz zCoC-i - ^ 2 | CoC-i + E [ CmCm-i + 2^Ti CtC?-i | » 

I m=l ~" J 

f zC-iCo - ^ 2 { C-iCo + E l^rrr^- iC + 2ZTi c ^+ 1 

I, m=l " 1 



dz ^OCl-t 



1 



(D-34) 

where Cm or C m stands for j m or n m . 
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The following polynomial of 1/z plays an important role in the calculations: 



Rm,k — Z ( n mjk jm n k 
[(ro-fc-l)/2] 

- E ( 



r=0 



(m — k — 1 — r) ! r (m + i — r ) 

If ^ 

r! (m — k — 1 — 2r)l r (k + § + r 



for m> k and 



m,k 



~Rk,: 



m— k~ 1 — 2r 

(D-35) 



(D-36) 



for m < k. By construction, a recursion formula similar to that satisfied by the 
spherical Bessel functions holds: 



R +R - 2fc + 1 P 



(D-37) 



Note that the indices m and fc can be negative as well. As examples, we write down 
the explicit forms for some special cases: 



Rk,k — 0, Rk,k+1 — 1, Rk,k-i - 

D.2.3. Source terms 

The source term can be rewritten as 



-1, 



Rk,k- 



2k -1 



(D-38) 



w\ n) =z 2 



+ 



' d 2 2d f 

dz 2 z dz V 

dz z 



£(£ + 1) 



Jn-l) 



2zj- + 1 
1 *(*»-!) 



(D-39) 



The contribution to Q from the first term is given by z _ C . So we focus on 
the second and third terms. Note that the operators of the second and third terms 
have opposite parities; the second term is odd while the third term is even under the 
transformation z —* — z. To perform the integration of these terms, we introduce the 
concept of the standard form of source terms as follows. 

First consider the first iteration, n = 1. Since = je and since we only need 

to calculate the real part of ^\ we only need to consider the second term. Using 
the recursion relations ( p-28j ) and (D-32), it may be rewritten in the form, 

otazjt + Poh-i + Mi+i ■ (D-40) 

Then using the integral formulas (D-34), the integrals F^^[zjg\ and i ? e </[j£±i] are 
readily evaluated to give 



Ft,i [zj t 
F tl [jk] 



D 



i-i 



R wo + E ( h + 

m=l 



rn 



1 



(£-k){£ + k + l) 



m + 1 
j k (k = £±l). 



Re,mjr, 



(D-41) 
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The real part of is expressed in terms of these functions, while the imaginary 
part is (]nz)j£ as given by Eq. (D-9). The result is Eq. (3-35) with q = 0. 

At the second iteration, n = 2, we insert the real part of Q to the second term 
in Eq. (D-39) and the imaginary part i(lnz)ji to the third term, to evaluate the real 
part of Q . Let us focus on the contribution of the terms of the form Re tm j m in ffl 
for the moment. Since Rt m are polynomials in 1/z, we cannot apply the integral 
formulas (D-34) directly. So, by using the recursion relation ( p-28| ) we get rid of the 
inverse powers of z. Then we find the corresponding source term may be expressed 
in the form, 

z(a-ji-i + a+ji+i) + Pmje+2m • (D-42) 

Similarly, at the third iteration, n = 3, the terms in f^p having the form z k j n (k < 0) 
will give rise to the source term which can be written in the form, 

z(a je + a_j'_£_2 + a.+3-l) + Y /W*+2m-i , (D-43) 



which is a generalization of Eq. ( D-40 ). Because the operator of the second term in 
Eq. (D-39) has the odd parity, the source terms for n = 2 and 3 take different forms. 
We call Eqs. (D-42) and ( D-43| ) the standard forms. For convenience we call the 



former the even standard form and the latter the odd standard form. Now turning 
to the terms with or (In z)j m , since they satisfy the same recursion relation as 
j m do, the same idea can be extended to them in a natural sense. The standard 
form for them is then defined by Eqs. ( D-42j ) and ( D-43| ) with j m replaced by or 



(\nz)j m . Note that D™ J at n = 2 plays an analogous role of jj> at n = 1. Hence the 
odd standard form of appears at n = 2. On the other hand, since the second 
and third terms in Eq. (D-39) have opposite parities, the parity of the standard form 
of (In z)j m is equal to that of j m . 

To summarize, the source term at the second iteration consists of the standard 
forms of 

j m : even, (\nz)j m : even, £>™ j : odd. (D-44) 
The integration of these terms can be done by using the formulas given in subsection 

(2) 

D.3 below. The resulting Q are given by Eqs. (3-40) for £ = 2, 3 and Eq. (4-16) for 
I = 4. Then we find there appear new types of the source term at the third iteration, 
which are 

zD™*, D™{, z(lnz) 2 j e , (In z) 2 j e±u z(\nz)D^£ 1 , (D-45) 



in addition to the opposite parity terms of Eq. ( D-44 ), 

j m : odd, (In z)j m : odd, D% : even. (D-46) 
D.3. Reduction of integrals 

In this subsection, we reduce the expressions [source terms] to those written 
in terms of F) J m . For this purpose we need to evaluate integrals such as 

F k/ [zDf] := n k J dz zj t D? J - j k J dz zn e Df. (D-47) 
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As an example let us show how this is evaluated. Using the basic integral formulas 
(D-34), we integrate it by part as 



F Kl [zDf\= n k J dz zj Df - j k J dz zn D% J 



3k 



-n k 



e-i 

noh + E 

m=l 



1 1 

+ 



m m + 1 
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-n k 



< 3o3o + 2^ — + TT 
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z 2 \J £/l 1 
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n m n m + -n e ne 
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-y) ^on +E (- + ^ 

/ k m=l 

/ ^ (- y ) {iono + E + ^1 ) ^ + ^ } ^0 



+Jfc 
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. . ^ / 1 1 

3o3o +Z^\ — + TT 

~i\m m + 1 

{n n + E ( ^ + 



m=l 



m m+1 



' x V m + 1 



3m3m + j3t3i > zn D 



n m n m + -n e ne > zj D 



n m3m 



=Dt nJ + \ UokD^ + E Q + J-) R mk D^ + - £ RekD 
+ \n k J zdz ( E d + T^y) Rmo3m + ]R^ Di 

~\jk I Zdz ( E + J^j) RraOUm + ^0^) itf 

\m=l / 

=or + \ + £ i (I + ^) w + 

k m=l x 7 



n J 



nJ 

i 



+ ~H k 
2 



E 



+ 



' V m m + 1 



Rmojm + -jRmje 



(D-48) 



The reduction of the term in the last expression is done similarly. In the following, 
we give formulas for each type of the source terms. 
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D.3.1. j m -terms 

The source terms have the standard form ( p-42 ) or ( D-43| ). 
Using Eqs. (D-34), their integrals are evaluated as 

F t,i [Cm] = -To UT-, 7~^C m (m t- £, 



F U [zQ] = Df 
Fe,e [zji-i] : 



(£-m)(£ + m + 1) 



1 

2 

DT 



Re,oCo + 5^ ( — I rrl R e,mCn 

m=l 



f-1 



4m 



-Rifln^ + E 4 m 2 _ ]_ Re,m3m- 



m=l 
-1 



4m 



[zn^_i] = Df - \ Rifijo + 4m 2 _ x R * 



771=1 



4m 



F?^ — —Df - -j Rg -i(o + E - — g — - -R^,m-lC. 



m=l 



(D-49) 
(D-50) 
(D-51) 
(D-52) 
(D-53) 



where C represents j or n. Note also that rather general formulas, 

1 

Fk,e[(m) = ■ 



. ,s {Rk,i(m+1 + Rk,£-l(m} 

m) {t + m + 1) 

" k,£ -( m (m^e, - 



m z 



1): 



if l ~ X / 1 1 \ 

K 2 '.. \m m + lj 



.Cm + -^RklCl 



(D-54) 
17-55) 



hold. 

D.3.2. (In z)j m -terms 

The source terms are in the form (|D • 4C| ) or ( |D • 42| ) with j m replaced by (In z^)j m . 
First we give the general formula for Fgi[(lnz)j m ]. Using the first formula in 
Eqs. (D-34), we obtain 



F^i[{\uz)j m ] 



(£-m)(£ + m + 1) 



In z 



1 



+ - 



1 



m 



-Ft, 



2z 



(£ + m+l) 



+ m + 1 



Jm+1 ~\~ Jr, 



(D-56) 



for m ^ £, — £ — 1. 

Next we consider the remaining term of the odd parity. With the aid of the 
second formula in Eqs. (D-34), after the integration by part, we have 

\ z { / 1 1 \ 1 

F ei [z(lnz)je] = F eo [z(lnz)j ] — < Re,ojo + Y] ( 1 — r ) Re,mjm > 

2 — ' \m m + 1 / 

k m=l x ' ) 



\m=l 



\m=l 
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~\ E {- + -4t) I 2 ( E r) + 7 \ n l ^,«J^D-57) 

4^Vm m + iy^ ^^fcy m(m + 1) j 
To evaluate the first term, we use the following trick. Note that 

Bjnj = [ dz zj (n Bjj - j B nJ ) 
Jo 

f z f z dz 

= dz zriQ (joBjj + n Q B nJ ) - / — B nJ 
Jo Jo z 

= B njJ - (In z)B nJ + fdzQnz)^-, (D-58) 
jo dz 

where we assumed J do not end with n, i.e., J / n,jn, ■ ■ ■. In the same way, 

rz dB j 

Bjjj = -B nnJ + (In z)Bj j- / dz(lnz) — 3 —. (D-59) 

JO "-2 



Thus we obtain 

'o 

dJ3< 



y ob(lnz)-^ = B jnJ - B njJ + (In z)B nJ , 
rz cm . 

/ dz(lnz)—^ = -Bjjj-B nnJ + Qnz)Bjj. (D-60) 

JO "-2 

Using Eqs. (D-60), we find 

F efi [zQnz)j ] = - DP + Qnz)D?. (D-61) 

As for the remaining terms of the even parity, we have 
F t)t [zQnz)j t+1 ] = -D™> + Df - \nzDf 

^ 4:777, 

-Re,-iQnz)jo ~ lnz z2 4m2 _ i R e,m-ljm 

^ <4?77 1 

+F t - 1 [zj ] + J2 4m 2 _ l F l,m-l[zjm] + ^—F e/ [zji +1 ], 
m=l 

FtM^)jt-i] = - D T j ~ On*W + \(^?k 

4m 



-R e _ 1 (]nz)j + In z ^ 4m 2 _ i R e,mjm-i 

m=l 

4m 1 
-F ti0 [zno] + ^^mfe-il + ^-Y^^NIA-]62) 

m=l 

D.3.3. (lnz) 2 j m -terms 

The source terms we have to evaluate are z{\nz) 2 ji and (\nz) 2 je±\. 
Using the first formula in Eqs. (D-34), we obtain 



F e/ [(lnzf Jm ]= {i _ m){e + m + 1 y 



/, n2 21nz 
-(lnz) 2 + 



£ + m + l (e + m + lf 
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(In z) (-2zj m+1 + (£ + m + l)j m ) + ' 



£ + m + 



T jj-,63) 



for m ^ £, —t—\. 

Also, with the aid of the second formula in Eqs. (D-34), we find 



F e/ [z(]nz) 2 j e ] =F efi [z(lnz) 2 j ] - I R e , j + ^ - + — — - ) R e , m j r 



2 i! ' ^— (. V. /// m + 1 

k m=l v 

€-1 



/ 1 1 \ 

+fi )0 [^(lnz)jo] + S — + — rr f ^ i z ( lnz )jm] 

m=l 

+-^[z(lnz)^]. (D-64) 

The evaluation of the first term in the above equation is done as follows. Using 
a technique similar to the one used to derive Eqs. (D-60), we obtain 

J dz (In z) 2 ^ = 2 [— B nnn j - Bjj n j + B n jjj — Bj n jj] 

+2 (In*) [B jnJ - B njJ ] + (lnz) 2 B nJ , 

J dz (lnz) 2 —j^- = 2 [-BjnnJ + B n j n j + Bjjjj + B nn jj] 

-2 (In z) [B nnJ + Bjjj] + (In zf Bjj. (D-65) 
Using these, we can rewrite the first term as 

F efi [z(lnz) 2 j ] = 2 (-£>™ nj ' + Dj jnj + + D} n ») 

-2(\nz) (D m + + (In z) 2 Df . (D-66) 

We need one more formula to evaluate (D-64): 

i^ m [z(ln2!)j m ] =F^ [2;(lnz)jo] - -^Rt, m jm 

In z I v-^ / 1 1 



fc=l 

m l\ 1 f / m 1 

1 I -^n -i 1 „ / X > 1 



4£(h^)H£j)^}w 7 ) 



D.3.4. _D£ J -terms 

The source terms are the even and odd standard forms of Djfy, and zD™ nj and 

n nnj 

u t±i- 

Necessary formulas are 

m=l 
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E 

vm=l 



-Rm,0jm-1 + r,„ , \ Je 



Am 2 — 1 



F^[zD^ 1 ]=Df J -R^Di J -Y, 
+H/ 



Am 



2t + r 



rn=l 



Am? - 1 



m J 



E4m . o . 

4m 2 _ 1 ^m-l,0Jm + 2£Z^^ 
m=2 y 
D nJ 



m)(£ + m + 1) 



.(£-m)(£ + m + l) 



(-Rm+1,0^ + ^m,0^-l) - ^ m ' J> 

t — m 



FtAzDf] 



F eA D e+i\ 



■ e-i 



E 



+ 



' V m m + 1 



£>#i + Hi 



nJ i 
£-ll 



2(e + iy 

■\nj 



2£ 



z 
2£ 



Rmfijm + -jRlfljl 



(Ro,e+2je + Ro,e+i3e-i) + Ro,e+de 

(D-68) 



2(£ + l) 

(Ro,dt + - Ro,e-de 



In order to evaluate the terms involving the integral operator Hf , we recast its 



argument into the form, 

z(a- 2 j-2 + a-ij-i + a j + ftiji) + ^ /3„j„. 

n^0,-l 

Then all the necessary terms can be easily evaluated as 

Hj [jrn 



(D-69) 



7 1 , ^ (--fym+ijo + Re, m no) + — ^j'o -1), 
m(m + 1) mz 



H\\zh\ = -DP + i2 <)0 no - R e ,ij , 
H 3 t [zj ]=D?, 
Hjlzj- 1 ] = -D^ 

Hj[zj-2] = -Df 1 - joRe-2 +j-iRe-i, 



(D-70) 



and 



! 



H™ J [zj-2] = (unnecessary for our present calculation), 

m(m + l) + ^.m-i^) 

-— 22 fc , m I>J J ' + 1 H{[zj m \. 
mz m(m + 1) 



H nj \i 



(D-71) 



The last term Hi[zj m ] can be reduced recursively to those given in Eq. (D-70). 
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D.3.5. z(]nz)D%-tenns 

What we need to evaluate are F t 
First we consider Fg g 



z{\nz)D n l i 1 ]. 
. It is evaluated as 



zQnz)Dfl 1 

zQnz)D?] - (lnz) iR^D^ + £ ^—R^ m _ lD 

K m=l 



m 



zu i+\ 



* 4m 1 
z(lnz) j E 4m 2_ 1 ^,oim-i + 2£Xi^+i,oi< 



where the first term is further expressed in terms of Df as 



Ft 



,-i 



s(ln«)£)J J 1 = -£>™ nj + nj - (In z)Z)f j 



The other unknown terms in Eq. (D-72) are also evaluated as 



~V,m— 1 



zDl j 



zi^ 



' m— 1 



fc=l 



4A; 2 - 1 



2m- 1 



(D-72) 
(D-73) 

Dl j \ 



4fc _ . i? mj0 



2m -I 



3 m— 1 



(m > 1), 



^ [z(lnz)io] = F AO [z(lnz)io] = -D™ - D\ nj + (lnz)D^ , 

Hj [(lnz)j m ] = (lnz) ( -. 1 ■ -r (-ifym+ijo + Re, m n ) + ^^jo) 
L m(m + 1) mz J 



2z 



m(m + 1 



;Jm+i 



2m + 1 

H — jTf (-fym+ijo - Re, m no) 



m 2 (m + 



,m . 

2 — JO- 

m z z 



(D-74) 



In the same way, we can evaluate Fg : g z(lnz)D"^ 1 j 



as 



e-i 



z(lnz)Df 1 ] = Fg fi [z(lnz)^] - (lnz) J Rg, D n _{ + ^J^^-i 

I m=l 



m=l 



zu i-\ 



I ipj 4 m i 

(lnz) ho + E 4m 2 _ 1 -Rm-l,0jm + ^^Rt-l,ojt 

(D-75) 

where the unknown terms in the right hand side are evaluated as 
Fgfi [z(lnz)D^l = -D? ji - D J g njj + {In z)D^ jj , 
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m—1 



Ft 



zD nj 



1,0 



m—1 



D?> - 



fc=i 

4fe 



4A- 



Rr 



2m -I 



m—1 



fc=i 



4/fc 2 - 1 



D . . Rm-1,0 . 

-Kfc-l.Ojfc + 7, 7 Jr. 

2m — 1 



(m > 1) 
(D-76) 



D.4. The asymptotic behavior 

Using the results of the preceding subsection, we obtain the analytic expression 

(3) 

for Q . The real part of it is given in Eq. (4-14) for 1 = 2 and in Eq. (4-15) for 
£ = 3, while the imaginary part is determined by Eq. (D-9). To obtain A™ c to 0(e 3 ), 
we then see that all what we need to evaluate are the asymptotic behaviors of D™ 1 , 
D™ nj , D™ nn:) and F£ t o[z (In z)jo]. Although the last of these can be expressed in terms 
of Df as given by Eq. @>6j]) , we find it is easier to evaluate the integral directly as 
it is. 

Here in order to evaluate the asymptotic behaviors of these functions, we first 
give necessary basic formulas. Then we evaluate the asymptotic behavior of all the 
necessary Bj and Fi i0 [z(lnz)jo]. In this subsection x = 2z. 

D.4.1. Basic asymptotic formulas 

First, we give the most basic formulas: 



S ■ 
C 



dy 



dy 



dy 

dy 
smy 

y 

cos y 



smy 

y 

cos y 



7T 

2' 



In x — 7, 



In y 



7T 



r7= 



y 

dy (my) 



Iny 



dy 



y 

cos y 



1 



II 



(\nyf 



IT 2 71-3 

2 7 + 24' 

--(In x) 3 - -L + TL. - -C(3), (D-77) 
3 V ' 3 12 3 sy 1 v ; 



where C,(z) is the Riemann zeta function and ("(3) = 1.202 

There appear several expressions to be estimated that diverge if we evaluate 
them term by term. But they give finite results when combined together. They are 

^2 



dy 

o y + l 

x dy 

o y + 

dy 



{Hv + 1) 



7T- 



Iny} -> — , 



o y + l 
y 



T {(ln( y + l)) 2 -(lny) 2 } 



{]n(y + 2) -Info + !)}-> — 



0, 

7T_ 2 

12' 



f ~~ ~l7 + 2 )) 2 - ( ln ^ + 1 )) 2 } ~+ 2 ^ 3 ' V2) + In 2 - ^(In2) 3 . 

Jo y + 1 L J 6 3 
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Jo 



x dy 



o y + 2 



ln(y + l){ln(y + 2)-ln(y+l)} 



„ . 1„ „. 3 3 



-0(3, 1/2) - — In 2 + -(In 2f + -C(3), 



(D-78) 



where 4>(a,b) represents the modified zeta function. We mention that the 0(3, 1/2) 
terms are found to cancel out in the final expression for A™. 

There are several formulas which require multiple integrations. They are evalu- 
ated as 

/ — / — (cos(n - y) - cosu) -> — , 
Jo y Uy u J 6 



[*dy \[°°du 

/ — my / — (cos(u — y) — cos u) 
Jo V Uy u 



7T 2 7 
IP 



(D-79) 
(D-80) 



f x dy \[°°du 
/ — lny / — (cos(u + y)- 
Jo y Uy u 



cos-u 



0(3,1/2) + ^ In 2- I(ln2) 3 + y ; 



f dy 
Jo 2/ 



/ — hiy / — sin u - j/ 

Jo y LA « 

f°° du. . , 
/ — (sm(u + y) 

[ x dy \ [°° du . . . 
/ — / — sm(u-y)- 
Jo y Uy u 



smti 



12' 



smti 



smti 



7T 

24' 

-► 0. 

u 



(D-81) 
(D-82) 
(D-83) 
(D-84) 



These are obtained by changing the variable from u to u' = — — 1, performing the 

1/ 

^-integration first and using the formulas (D-78). 

Next we give several formulas that contain S and C. Let us recall the definition 
of S and C: 

n f x dy . f°° dy . ir 

S= — smy = - — smy + -, 

Jo y Jx y 2 



dy 



f x dy f c , 

C = / — (cos y — 1) = — / — cos y — In x — 7. 
Jo y Jx 



y Jx y 

By using the integration by part, we have 



(D-85) 
(D-86) 



rx 

/ dy lny 
Jo 



sin y cos y — 1 



y 



-s + — - — c ' 

= 7 L [ dy^^-lny- [ dy^^- (7 lny + (lny) : 

2 Jo y Jo y v 

rX y rOO ^y 

— I dy / — (cos(u — y) — cosn) 

Jo y Uy u 



cos y — 1 



(D-87) 
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whose asymptotic behavior is determined by Eqs. (D-77) and ( p-80 ). Similarly we 
have 



and 



dy lny 

o v y 

7T 



siny _ cos y - 1 



sin y / - cos y - 1 

lny + / dy 

2 Jo y Jo y 

x lny f f 00 dii 



(7 lny + (lny) 2 ) 



+ / dy- 
Jo 



y 



•i 



u 



- (cos(u + y) — cos u) 



(D-88) 



o V y y 



j sm f i n i \ I I y] cosy-1 

dy lny(lny + 7) ± - / dy lny 

io y 2 Jo y 

rx In y f f 00 dit , 



Z* 3 ' In y /*cx^ Qf-^ 

/ dy / — (sin(it ± y) — sin it) 

Jo y Uy u 



(D-89) 



whose asymptotic behaviors are determined by Eqs. (D-77), (D-81), ( p-82| ) and 
( jDjgj ). Then we obtain 



f ^ (s 2 + 
Vo y v 



^(lnx + 7 ) + 5 (In a: + 7 ) 3 - ~C(3), (D-90) 



f *L (s 2 - & 
Jo y v 



7T 2 1 Q 7T 2 T 4 

> T lnx--(lnx + 7 ) 3 + T I + 3C(3) 
-2U(3,l/2) + ^ln2-i(ln2) 3 j . 



(D-91) 



These two formulas are obtained by reducing them to the forms to which Eqs. (D-87) 
and (D-88) can be applied, respectively. 

Finally, we present two complicated formulas. The first is 



dy 
o y 



-2 sin ySC + (cos y - 1) (S 2 - C 2 



>2 ^(3,l/2) + ^ln2-l(ln2) 3 j 
-gC(3) - ^vr 2 (lnx + 7) + ^(Inx + 7 ) 3 



(D-92) 



where we have used the equalities 
sin(t + 2)x 



SC 



S 2 -C< 



I _l_ ln(t + l)dt + ^C-(lnx + 7 )5 + |(lnx + 7), (D-93) 

°° COS(t + 2)X ln(t + l)dt + ttS - ^ + 2(lnx + 7 )C + (In x + 7 ) 2 , 

(D-94) 
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and applied the formulas (D-77), (D-87) and the last one of Eqs. (D-78). The second 
is 

x dy 
io y 



J 2 : = f X JL (Vcos y - 1)SC + smy(S 2 - C 2 ) 
Jo V v 



7T 
'2 



(lnx + 7 ) 2 + ^), (D-95) 



which is obtained in the same way by applying the formulas (D-77), ( p-84 ) and 
(D-89). 

D.4.2. The asymptotic behavior of Bj 

As we have mentioned, what we have to evaluate are the asymptotic behaviors 
of D™ J , D™ nj , D™ nnj and Fz j o[z(lnz)jo]- Hence, recalling the definition of Df, we 
need to evaluate B n j, Bjj, B nn j, Bj n j, B nnn j and Bj nn j in addition to F^ [-z(ln z)jo\. 

The formulas for Bj with two indices are given by the first two equations of 
(D-77): 

B nj = ~S^~, Bjj = -±C^±Qnx + i). (D-96) 
As for F( i o[z(lnz)jo], its asymptotic behavior is directly evaluated as 

Fi [z(lnz)jo] = n £ [ dzz(ln z)j$ - j £ [ dzz(lnz)j n 
Jo Jo 

- ~ {{hi zf + ^ - ( 7 + In 2) 2 ^j n e - ~( 7 + In 2)j t -(D-97) 
The formulas for Bj with three indices are given by 

J J A Jo V y y 

-£(lnx + 7 ), (D-98) 

o 

where we have used Eqs. (D-77) and ( p-84 ), and 

1 f x /cosy + 1 siny 

Bnnj = -7 / dy\ C H ; 

A Jo V y y 

-» -- f— vr 2 - (lnx + 7 ) 2 V (D-99) 



.12 

where we have used Eqs. (D-77). 

For Bj with four indices, we have 

^ Jo y ^ Jo y 

= -\ [SB jnj - CB nnj ] + l -[ X dy (^LC - C -^±S ) S 
z 8 Jo \ y 

1 r d ^ ig+ cosy+i c \ c 

8 7o V y y / 
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C 2 

= -7; [SBjnj ~ CB nnj ] + ; I dy— 



1 

24 



y 



V(lnx + 7 )-i(lnx + 7 ) 3 -C(3) 



(D-100) 



where we have used Eqs. (D-90), (D-91) and (D-92), and 



B 



nnnj 



L 



sin v „ 
dy -B, 



1 



y nn] + 2 
1 



dy 



cos y + 1 



B 



y 



dyBjnj + - [CBj n j + 5S„ n j] 
2 



1 



/" sin y cos y 
o — 



liaxB 



J n J 



V y 
1 



siny 



5 C + - 



dy 



c _cosy-l 5 



V y y 



1 1 1 r x SC 

+ - [C5 jnj + SB nnj ] + ^h + - A J o dy— 



IT 

32 



-7T 2 + (lnX + 7) 2 



(D-101) 



where we have used Eqs. (D-89) and (D-95). We should note that Eq. (D-89) has 

, ■ 1 f x . SC 

been also used m evaluating the term - / dy . 

4 Jo y 

With these results, the asymptotic incoming amplitudes Af c to the required 
order are obtained, which are given in Eqs. (4-17) in the text. 



Appendix E 

5. 5PN formulas for Rf and {dE/dt)^ m in the Schwarzschild case 



In this Appendix, we show the post-Newtonian expansion of R m which are 
necessary to evaluate the 5.5PN gravitational wave luminosity in the case of a 
Schwarzschild black hole. Then we show each (£, m)-mode contribution to the total 
luminosity from a particle in circular orbit. 

For convenience, we give the formulas for cqujR" 1 , where we recover the index I 
on R m and c = {I - !)£(£ + !){£ + 2)- 6te. 



cqujR\ 



4z A 8i * 



6 



24948 



22 z 
105 ~ 

59 z 12 



2i 



23 z h 



35 Z + 1890 ' 945 



+ 



2i 



.9 



13 z 



10 



+ 



„13 



41580 

83 z 1A 



v 15 



12972960 2162160 



1945944000 277992000 



+ 



3i 



5 5 
9403 z 9 293 i 



3118500 594000 



lz 5 13 i a 109z 7 341 i , 

z 6 H 1 z s 

63 90 1890 22680 



10 38963 z 11 75529 i 
z 10 + —=— + 



v 12 



567567000 9081072000 
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4z 2 123317 z 4 231479 i 5 889954 z 6 454499 i 7 
1 "5" + 36750 + 110250 * ~ 1157625 ~ 2315250 * 

215321483 z 8 35106811 % 9 214 z 4 In 2 428 i 5 
+ 5501034000 + 5501034000 Z 525 1575 Z UZ 

1177 z 6 Inz 107 i 7 2461 z 8 Inz 107 i „ , . , 

+ 11025 + 3675 Z 396900 " 99225 "' 

. -66823 z 3 99851 i 4 504569 z 5 2488639 i 6 428 z 3 Inz 
1 12250 55125 * " 694575 ~ 3969000 Z + 525 

107 i 4 , 428 z 5 Inz 1391 i 

H z In z H 1 

350 6615 18900 

/ 471487 z 2 263i 3 214z 2 lnz\ , 
+ [ 220500 " 1260 Z 525 1 ' ' 

nlll , 4z 5 2i fi 2z 7 i 8 29 z 9 i 10 47 z 11 
c ujRf = + — z 6 z 8 + — — + z 10 - 




21 21 63 135 20790 4620 1621620 

z 12 + Ji z !L + '^ z u 



294840 64864800 29937600 

. -10z 4 53i = z 6 i 7 751z 8 1483i q 

+ z 5 H z H 1 z 9 

" 21 315 210 90 155925 1247400 

23 z 10 367 i 

■ z 



102960 10810800 

8z 3 i 4 40337 z 5 79099 i 6 12562147 z 7 4840537 i g 
+ '~2T + 14 Z + 46305 + 185220 Z 91683900 ~ 157172400 Z 

26 z 5 Inz 13 i R , 13 z 7 Inz 13 i . ; 

z In z H 1 z 8 In z e 2 

441 441 1323 5670 

/ -2z 2 182981 z 4 3753697 i 5 

+ ^~21 92610 5556600 Z 

65 z 4 Inz 689i 5l \ 3 

H 1 z 5 Inz e 3 , (E-2) 

441 13230 y ' V ; 

nlll , 2z 6 U 7 13z 8 8z 9 
" 3465 " 10395 " 

71 z 10 . i „ 37 : 2 , a , 



+ -: r « -: ; z 1 



540540 54054 16216200 4054050 
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4i fi 142 z 7 31 i s 929 z 9 o, ln 
135 51975 51975 2702700 96525 

5 958223891 z 6 239560304 i 7 
Z + 6051137400 + 3781960875 Z 

1571z e lnz 3142 i 



c ujR 



in 




(E-3) 



(E-4) 



19305 135135 135135 405405 
-32 z 1 163 i s 




z* e, (E-5) 



cqujR 1 ? = + z w . (E-6) 

7 225225 225225 v ' 

Next we show the contribution of each (£, m)-mode to the gravitational wave 
luminosity to 0(v n ) in the case of a circular orbit around a Schwarzschild black 
hole. We set 

/dE\_ 1 (dE^ 



^2ve,m, (E-7) 



where (dE/dt)-^ is the Newtonian quadrupole luminosity, Eq. (4-19). In the present 
case we have rj£ m = T]i- m - Hence we show only the modes m > 0. 

107 v 2 o 4784 v 4 428 ttv 5 

572,2 = 1 ^ h47lV + 



21 1323 21 

6 / 99210071 1712 7 16 vr 2 3424 ln 2 1712 InV 

+ V \ 1091475 105 + ~~ 3 105 105 / 

19136 ttv 7 



1323 



27956920577 183184 7 1712 tt 2 366368 ln 2 183184 Inu 



+ V [ 81265275 + 2205 63 + 2205 + 2205 



100 



+ v y 



( 396840284 tt 6848 7 vr 13696 vr In 2 6848 vr lnv 



1091475 105 105 105 



10 / 187037845924 8190208 7 76544 vr 2 
+ V \ 6257426175 138915 + 3969 
16380416 In 2 8190208 In v s 



138915 138915 

, , / -111827682308 vr 732736 7 vr 

+ v 11 1 — 

81265275 2205 

1465472 tt In 2 732736 tt In v\ 

+ ™ . (E-8) 



v 



2205 2205 

2 17 V 4 7TV 5 2215 v 6 177TV 7 



V2,1 ~ 36 ~ 504 + 18 254016 252 

8 / 15707221 107 7 vr 2 107 In 2 107 1mA 
+ V ^26195400 ~ 945 + 27 945 945 J 

2215 7r 
~ 127008 

w ( 6435768121 1819 7 17 vr 2 1819 In 2 1819 1mA 
+ V \ 57210753600 + 13230 ~ 378 + 13230 + 13230 ) 

u / 15707221 tt 21477T 214 vr In 2 214vr lnv\ 

+ V V 13097700 945 945 945 ) ' ( ' ' 

_ 1215 v 2 1215 v 4 3645 tt 243729 v 6 3645 ir v 7 

m ' 3 ~ 896 112 + 448 + 9856 56 

o / 25037019729 47385 7 3645 vr 2 47385 In 2 

+ v 8 H 

^ 125565440 1568 224 1568 

47385 In 3 47385 Inv' 



+ 



1568 1568 
7311877TW 9 



4928 



in / 2074855555161 47385 7 3645 vr 2 47385 In 2 

+ u 1 1 

^ 1381219840 196 28 196 

47385 In 3 47385 \nv 
+ 196 + 196 

n / 75111059187 tt 14215577T 142155 vr In 2 



+ v 



■ 



62782720 784 784 

142155 vr In 3 142155 vr In v s 



784 784 
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Appendix F 

4PN luminosity in terms of the orbital frequency 



Here we present the (£, m)-mode contributions to the gravitational wave luminos- 
ity for a circular orbit on the equatorial plane around a Kerr black hole. In stead of 
v = (M/ro) 1 / 2 , the formulas are expressed in terms of the parameter x = (MI?^) 1 / 3 , 
where f2 v is the orbital angular frequency, which is more relevant in the actual analy- 
sis of observed gravitational wave signals. We express the partial mode contributions 
as 
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Appendix G 

4PN luminosity for a slightly eccentric orbit in the Schwarzschild 



case 



We define the partial mode contribution t]t^ n 



as 



dE\ _ 1 (dE\ 

~dt/~ 2 ¥ N .^ % ' m ' n ' 

1N £,rn,n 



(G-l) 



where (dE/dt)^ is the Newtonian quadrupole luminosity, Eq. ( (4-19|) , and r\i_ m _ n = 
Ve,m,n because of the symmetry in the Teukolsky equation. Up to 0(t> 8 ), do 
not appear irrespective of £, and other modes are as follows. 
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Appendix H 

Asymptotic amplitudes and R up 

In this Appendix, we show the asymptotic amplitudes B mc , B tra,ns and C tia,ns 
and the post-Newtonian expansion of i? up which are used in section |l2| to evaluate 
the black hole absorption rate to 0(v 13 ) relative to the quadrupole energy flux at 
infinity in section, 
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Appendix I 

Energy absorption by a Kerr black hole 

Here we give the (£, m)-components of the energy absorption rate to 0(v 8 ) 
beyond the lowest order for the Kerr black hole, that is 0(v 13 ) relative to the 
quadrupole luminosity at infinity. 
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Appendix J 

(dE/dt)n in terms of the orbital frequency 

In this Appendix, we describe the absorption rate (dE/dt)u by a Kerr black hole 
in terms of x = (M^) 1 / 3 . Using the relation, 

V = x(l + ^x 3 + ^ 2 x 6 + 0(x 9 )), (J-l) 

we have 
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